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ABSTRACT 


A l/12th-degree, pan-Arctic ice-ocean numerical model is used to better understand the 
circulation and exchanges in the Bering Sea. Understanding the physical oceanography 
of the Bering Sea is significant for the U.S. Navy due to the expected increase in ship 
traffic and exploration of natural resources that will likely coincide with the ongoing 
retreat of sea ice in the Western Arctic. This model represents a large step forward in the 
ability to simulate the mesoscale eddies and meanders in the Alaskan Stream and the 
deep Bering Sea basin, which are shown to exert a strong control on the flow into and out 
of the western Aleutian Island passes. Model results show that upwelling of deep Bering 
Sea water, which is the primary source of nutrients for important ecosystems of the 
Bering, Chukchi, and Beaufort seas, is enhanced by the presence of cyclonic eddies in the 
vicinity of canyons along the slope. High values of eddy kinetic energy in Bering and 
Anadyr straits help explain the areas of high biological productivity located just 
downstream in the Chirikov Basin and north of Bering Strait. Model results show 
significant horizontal and vertical shear in the flow through Bering Strait, and indicate a 
need for more observations of the flow structure on a continuous basis. 
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I. INTRODUCTION 

A. GEOGRAPHY OF THE BERING SEA 

The Bering Sea is a high-latitude, semi-enclosed, marginal sea, which links the 
North Pacific Ocean and the Arctic Ocean (Figure 1). To the south, the Bering Sea is 
separated from the North Pacific Ocean by the Aleutian Island Arc, a chain of 
volcanically-formed small islands and passages. The Bering Sea is composed of 
approximately equal areas of continental shelf and deep basin. The deep basin makes up 
the southwestern part of the Bering Sea and has a maximum depth of -4,000 m. Two 
submarine features (Shirshov Ridge and Bowers Ridge) create three sub-basins within the 
deep Bering Sea: Aleutian (the largest), Kamchatka, and Bowers basins. The Bering Sea 
slope is -1300 km in length and runs from the easternmost Aleutian Islands to just south 
of Cape Navarin, Siberia. 

The northeastern portion of the Bering Sea is continental shelf, which is extensive 
in the east (>500 km) and narrows to the west (<100 km). Most of the Bering Sea shelf is 
<100 m deep, and the shelf break occurs at -170m deep (Okkonen 1996). Three straits 
are located on the northern Bering Sea shelf: Anadyr Strait (between Siberia and St. 
Lawrence Island), Shpanberg Strait (between St. Lawrence Island and Alaska), and 
Bering Strait (between Siberia and Alaska). The Bering Strait is the northernmost border 
of the Bering Sea and is the only Pacific connection to the Arctic Ocean. 
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Figure 1. Model bathymetry (m) and geography of the Bering Sea. 
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B. MAJOR FEATURES OF THE BERING SEA 

The Bering Sea, like many subpolar seas, is characterized by seasonal extremes in 
solar radiation, meteorological forcing, and ice cover (Stabeno et al. 1999). The primary 
storm track in the North Pacific is along the Aleutian Island Arc, which forms the 
statistical feature known as the Aleutian Low (Overland 1981; Overland and Pease 1982). 
During summer, the Aleutian Low is usually weak and weather is mild, due to the long 
periods of daylight and high insolation. Significant changes occur in atmospheric 
pressure fields during winter, with a strengthened Siberian High and an intensified 
Aleutian Low. This results in strong, very cold winds from the northeast and an increase 
in the frequency and intensity of storms. 

Portions of the shelf and coastal region experience seasonal sea ice cover for 
approximately nine months per year (October-June), with ice reaching the shelf break in 
late winter of some years, although satellite observations indicate a decrease in the extent 
of sea ice and a lengthening of the ice-free period since 1979 (Frey et al. 2011). 
Maximum ice extent tends to occur in March, with little to no ice present from July- 
September. Active polynya regions include the St. Lawrence Island Polynya (just south 
of the island), as well as coastal areas just south of the Siberian Peninsula and Seward 
Peninsula (Schumacher et al. 1983). Intense ice fonnation in these areas leads to 
salinization of the water column and the erosion of the thermocline (Schumacher et al. 
1983). Most of the sea ice in the Bering Sea is formed in situ, with very little moving 
southward through Bering Strait (Tabata 1974; Muench and Ahlnas 1976). Large-scale 
ice motion (predominately southward) is driven primarily by regional wind forcing, 
except in local areas where currents may be strong enough to oppose the wind (Muench 
and Ahlnas 1976). Pease (1980) described the sea ice process in the Bering Sea as a 
“conveyor belt” system, where ice growth occurs primarily in the north, ice is advected 
southward due to wind stress, and ice melt occurs at the thermodynamic limit. 
Subsequently, as meltwater cools the upper layer, the sea ice extent is allowed to advance 
farther south. 
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Major circulation features include the Alaskan Stream, a westward-flowing 
current located just south of the Aleutian Island Arc. The stream provides the inflow into 
the Bering Sea via the Aleutian Island Passes (Favorite 1974). The Bering Slope Current 
flows along the northeastern boundary of the Aleutian Basin above the continental slope 
(Kinder et al. 1975), although more recent observations (e.g., Okkonen 1993) and model 
results presented herein indicate that this current may be more a system of eddies rather 
than a continuous feature. The Anadyr Current is located in the northwestern Bering Sea 
and flows along the western boundary, generally toward Bering Strait. The Alaska 
Coastal Current (ACC) is a narrow (order of ten km; Woodgate and Aagaard 2005), 
buoyancy-driven current which originates south of the Aleutian Islands, flows through 
Unimak pass and other eastern passes, and continues northward along the Alaskan coast. 
Observational results indicate that the presence of the ACC may be seasonal at Bering 
Strait (Paquette and Bourke 1974; Ahlnas and Garrison 1984; Woodgate and Aagaard 
2005). 


C. SIGNIFICANCE OF THE BERING SEA 

The Bering Sea shelf is a highly productive region of the world ocean (e.g., 
Grebmeier et al. 1988 and references therein). Many species of fish, shellfish, seabirds, 
and marine mammals are sustained by large populations of the lower trophic levels (e.g., 
phytoplankton and zooplankton). The northern Bering Sea supports a historically 
abundant benthic population of large clams and polychaetes. These benthic organisms 
support populations of Pacific walrus, diving sea ducks, and grey whales. However, 
recent declines in these benthic populations may be affecting the higher trophic 
organisms (Grebmeier and Dunton 2000; Grebmeier et al. 2006), some of which are 
believed to have critically small populations. 

Almost half of the total U.S. annual catch of commercial fish and shellfish comes 
from Alaskan waters (Macklin et al. 2002; Sigler et al. 2010), primarily in the southeast 
Bering Sea. Important commercial species include pollock, salmon, halibut and crab. 
Because this region is an important national and international fishery, much fisheries 

research has been conducted in some parts of the Bering Sea (e.g., Busby et al. 2005; 
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Ciannelli and Bailey 2005; Mueter et al., 2006). For example, Sigler et al. (2010) note 
that walleye pollock experienced low recruitment during the “warm” years of 2001-2005, 
possibly due to very low abundance of their main prey (zooplankton). The pollock 
fishery, which had a wholesale value of $1.4 billion in 2008 in the U.S., was 
subsequently affected by decreased commercial catches (Sigler et al. 2010). Mechanisms 
for the climate-related changes are still not completely clear and are currently being 
investigated. In general, there is a need for more understanding of the physical processes 
that control environmental conditions for these organisms. This knowledge is important 
not only for scientific understanding, but also for people who rely directly on the Bering 
Sea for food and employment. Limited observations exist for the Bering Sea, as 
compared to temperate waters, and even less work has been done in other marine science 
fields, particularly in numerical modeling. Model results presented herein will advance 
understanding of important physical processes in the Bering Sea. 
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II. MODEL DESCRIPTION 


Early models of the Bering Sea were two-dimensional and barotropic (e.g., 
Overland and Roach 1987; Spaulding et al. 1987). These models were able to represent 
important features of the northern Bering Sea circulation. Overland and Roach (1987) 
showed that the large-scale (Pacific-Arctic) pressure gradient was responsible for driving 
the background northward flow through Bering Strait. However, their ability to represent 
baroclinic flow and shelf-basin exchange was restricted. More recent work in the Bering 
Sea includes that of Overland et al. (1994). They used a primitive equation, three-layer 
hydrodynamic model at 1/8° horizontal resolution. However, the domain did not include 
the shelf and slope, but was instead limited to the Bering Sea basin only. Mizobata et al. 
(2006) used a numerical model based on the Princeton Ocean Model (Blumberg and 
Mellor 1987) to simulate the eddy field in the vicinity of Zhemchug Canyon. While this 
model has a high horizontal resolution (5 km), the domain covers only a fraction of the 
Bering Sea slope and results may be significantly influenced by the prescribed lateral 
boundary forcing. The following analyses utilize results from a numerical model with a 
pan-Arctic domain that includes the entire deep basin, slope and shelf of the Bering Sea, 
as well as the North Pacific Ocean southward to ~30°S. 

This model was developed by the Naval Postgraduate School Arctic Modeling 
Effort (NAME; http://www.oc.nps.edu/NAME/name.html). The model is configured on 
a 1/12° rotated spherical coordinate grid (Arakawa B grid; Mesinger and Arakawa 1976). 
The coupled sea ice-ocean model has 45 z-coordinate vertical depth layers with eight 
levels in the upper 50 m and fifteen levels in the upper 200 m (Figure 2). The high 
vertical resolution, especially in the upper water column, allows for more realistic 
representation of the shallow Arctic and sub-Arctic shelves. The ~9-km horizontal model 
grid permits calculation of flow through the narrow straits of the northern and southern 
Bering Sea and permits eddies with diameters as small as 36 km, however, the smallest 
eddies are likely not resolved (Maslowski et al. 2008b). According to Chelton et al. 
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(1998), the typical Rossby radius of deformation range is -12-20 km in the Bering Sea. 
Thus, fully eddy-resolving models of the Bering Sea require a grid cell size of order few 
kilometers. 

The model domain contains the sub-Arctic North Pacific (including the Sea of 
Japan and the Sea of Okhotsk) and North Atlantic Oceans, the Arctic Ocean, the 
Canadian Arctic Archipelago and the Nordic Seas (Figure 2). The region of interest, the 
Bering Sea, is therefore far away from the artificially closed lateral boundaries in the 
North Pacific at 30°N, greatly reducing the potential effect of boundary conditions. In an 
effort to balance the net flow of Pacific Ocean water into the Arctic Ocean, a U-shaped 
500-m deep, 162-km (18 grid point) wide channel was created through North America 
connecting the Atlantic Ocean to the Pacific Ocean. A westward wind forcing of 1.75 
dyne cm" is prescribed along the channel. Flow through the Bering Strait and the channel 
is predicted in the model at each time step. Model bathymetry is derived from two 
primary sources: ET0P05 at five-minute resolution for the region south of 64°N and 
International Bathymetric Chart of the Arctic Ocean (Jakobsson et al. 2000) at 2.5-km 
resolution for the region north of 64°N. 

The regional ocean model adapts the Los Alamos National Laboratory Parallel 
Ocean Program model (POP; Dukowicz and Smith 1994). The combination of the free 
surface and high resolution allows the use of unsmoothed, very realistic bathymetry. 
This is important for representing steep bathymetry, such as continental slopes. Yukon 
(and other Arctic) river runoff is included in the model as a virtual freshwater flux at the 
river mouth using the observed daily-averaged annual cycle of each river’s discharge. In 
the Gulf of Alaska the freshwater flux from runoff (Royer 1981) is introduced by 
restoring the surface ocean level (of 5 m) to climatological (Polar Science Center 
Hydrographic Climatology; PHC; Steele et al., 2001) monthly mean temperature and 
salinity values over a monthly time scale (as a correction tenn to the explicitly calculated 
fluxes between the ocean and overlying atmosphere or sea ice). This approach was used 
in part due to the lack of realistic discharge observations for the many small rivers, which 
empty into the Gulf of Alaska, as well as a lack of domain-wide and time-varying 
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precipitation/evaporation data. A 4°-wide band of ocean points along the domain 
boundary is restored to annual average PHC temperature and salinity climatology on a 
ten-day timescale. 

The sea ice model uses viscous plastic ice rheology (Zhang et al. 1999; Zhang and 
Hibler 1997) and the zero-layer approximation for heat conduction through ice after 
Semtner (1976). The surface heat budget follows Parkinson and Washington (1979). 
The sea ice model was initialized on January 1, 1979 with a 2-m-thick slab of ice placed 
over all grid cells with sea surface temperatures < 0.0°C. 

The ocean model was initialized with climatological, three-dimensional 
temperature and salinity fields (PHC) and integrated for 48 years in a spin-up mode. 
During the spin-up, daily averaged annual climatological atmospheric forcing derived 
from 1979 to 1993 reanalysis from the European Centre for Medium-Range Weather 
Forecasts (ECMWF) was used for 27 years. The spin-up was continued using repeated 
1979 ECMWF annual cycle for six years and then 1979-1981 interannual fields for the 
last 15 years of spinup. This approach is especially important in establishing realistic 
ocean circulation representative of the time period at the beginning of the actual 
interannual integration. At the same time, the spin-up procedure was designed to force 
the model into a quasi-equilibrium state that is minimally sensitive to the specific initial 
conditions. The ocean time step is 480 s and the ice model time step is 48 min. The 
ECMWF atmospheric forcing fields include: 10-m east-west and north-south (u and v) 
wind velocity components, surface pressure, temperature and dew point, and incoming 
longwave and shortwave radiation. The final run with realistic daily averaged ECMWF 
interannual forcing starts in 1979 and continues through 2004. Results from this 
integration (26 years) are used for the following analyses. 
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Figure 2. 1/12° ice-ocean model domain and bathymetry (m). 
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III. ORGANIZATION OF THE DISSERTATION 


The following four chapters are organized geographically, beginning at the 
southern boundary of the Bering Sea, with Chapter IV on the communication between the 
deep Bering Sea and the Western Subarctic Gyre, as well as the circulation and 
mesoscale activity within this region. Chapter V is focused on the shelf-basin exchange 
and outer shelf dynamics. In Chapter VI, the circulation and energetics are examined 
across the northern Bering Sea shelf. The flow and property exchange through the 
northern boundary of the Bering Sea (Bering Strait) are intercompared from observations 
and multiple model results in Chapter VII. Chapter VIII shows the linkages of the Bering 
Sea modeled circulation with large-scale weather indices, as well as the physical impacts 
on biological processes, and the downstream effects of Bering Sea processes. The overall 
conclusions are discussed in Chapter IX. Finally, some suggestions for future research 
are presented in Chapter X. 
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IV. THE DEEP BERING SEA AND WESTERN SUBARCTIC 

GYRE COMMUNICATION 

A. INTRODUCTION 

The deep Bering Sea basin is bordered by the Aleutian Island Arc to the south, by 
the Bering Sea slope to the northeast, and by the Kamchatka Peninsula to the northwest. 
Depths greater than 3000 m are found throughout the basin, except in the vicinity of two 
large submarine features: Bowers Ridge (north of Amchitka Pass) and Shirshov Ridge, 
which extends southward from Kamchatka Peninsula. The model bathymetry of the deep 
Bering Sea, along with a series of cross sections to be discussed in the following section 
is shown in Figure 3. The general circulation of the Bering Sea basin is typically 
described as cyclonic in the long-tenn mean. However, transport within the gyre can 
vary by more than 50% (Stabeno et al. 1999). Causes associated with this variation have 
been identified as either changes in the Alaskan Stream inflow (Overland et al. 1994) 
and/or variability of the wind-driven transport within the basin (Bond et al. 1994). The 
low-frequency variations in the wind stress curl account for -1.5 Sv of the variability in 
the transport within the deep Bering Sea basin (Bond et al. 1994); however, the majority 
of the variability is due to variations in the curl on daily to monthly timescales. 

This chapter includes discussion of the general circulation within the Bering Sea 
basin and along the Aleutian Islands from the long-tenn model mean results. In addition, 
volume and property flux values across Aleutian Island Passes are shown and compared 
to available observations. Finally, the mesoscale meanders and eddies within the basin 
and along the Alaskan Stream are discussed and validated with limited observations. 
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Figure 3. Model bathymetry (m) of the deep Bering Sea. Black lines indicate the 

locations of cross-sections and white text indicates the names of these cross- 

sections. 


B. GENERAL CIRCULATION 


The Alaskan Stream flows primarily westward just south of the Aleutian Island 
Arc (Figure 4). This current is a major feature of the northern North Pacific, as indicated 
by the high values of total kinetic energy (TKE) in Figure 4. TKE is calculated as: 


TKE = 



( 1 ) 


with u and v representing the x- and y-components of velocity. TKE is useful for 
determining the relative speed of a current. The Alaskan Stream has a significant effect 
on the flow and property flux through the Aleutian Passes and, therefore, on the Bering 
Sea. The mean modeled volume transport is between 34 and 44 Sv, with intensification 
occurring downstream (Maslowski et al. 2008a). 
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Reference Vectors 
(cm/s) 

Figure 4. The 26-year mean (1979-2004) upper 100 m circulation (vectors) and total 
kinetic energy (shading). Magenta contour lines represent bathymetry (m). 

Every second vector is shown. The direction of the Alaskan Stream is 

westward. 

Communication between the Bering Sea and the North Pacific (primarily the 
Alaskan Stream) occurs through several passes along the Aleutian Islands. Inflow 
through these passes contributes to the eastward-flowing Aleutian North Slope Current 
(ANSC; Figure 4), which is located north of the Aleutian Islands. The ANSC is partially 
deflected by the bathymetry of Bowers Ridge and flow through Amchitka and Amukta 
passes. These are areas with relatively strong northward inflow, which tends to push the 
ANSC northward. The generally eastward flow of the ANSC begins to turn 
northwestward between 177-167°W, forming the Bering Slope Current (BSC). The BSC 
extends from the eastern Aleutian Islands along the shelf-break toward the coast of 
Russia. Long-term mean volume transport was estimated to be ~2-3.5 Sv (Clement 
Kinney et al. 2009), however the BSC is more a system of eddies, rather than a 
continuous current (Okkonen 1993; Clement Kinney et al. 2009). A major portion of the 
BSC separates from the slope around 180°, and gradually turns southwestward. The 
shallower portion is bathymetrically steered southwestward along the Siberian coast 
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within the 500-m isobath, fonning the Kamchatka Current. West of Shirshov Ridge, 
these two merge and flow primarily southward into the North Pacific Ocean. These 
general circulation features compare favorably with available observations published by 
Stabeno and Reed (1994), Verkhunov and Tkachenoko (1992), Reed et al. (1993), and 
Cokelet et al., (1996). 

C. VOLUME AND PROPERTY FLUX ACROSS THE ALEUTIAN ISLAND 

PASSES 

According to Stabeno et al. (1999), volume transport into the Bering Sea is highly 
variable on time scales of weeks to years. Inflow from the North Pacific into the Bering 
Sea occurs primarily through Near Strait, Buldir Pass, Amchitka Pass, and Amukta Pass 
(Stabeno and Reed 1993). Outflow from the Bering Sea into the North Pacific occurs 
largely through Kamchatka Pass (Stabeno and Reed 1993). Table 1 shows the model 
results of 26-year-mean (1979-2004) volume, heat, freshwater, and salt flux across each 
pass or strait shown in Figure 3. Heat flux is referenced to the freezing temperature 
(based on salinity) and freshwater flux is referenced to 33.8. Some passes were labeled 
alphabetically because they either did not have a recognized geographic place name or 
were artifacts of the 9-km grid cell spacing. The smallest islands were not resolved by 
the model bathymetry and at times created very shallow passes that are not present in 
reality. These alphabetically-labeled artificial passes had very low mean volume 
transport values, with most < +/- 0.006 Sv (the value at Strait C is -0.136 Sv). Model 
results have previously been compared with observations in Amuka Pass and Amchitka 
Pass (Maslowski et al. unpublished). Other passes in the central and eastern Aleutians 
also tend to have net northward flow, however they are much shallower and narrower, 
which means that they contribute little to the overall volume and property flux (Table 1). 
Therefore, this chapter will focus on the larger and deeper western passes including: 
Buldir Pass, Near Strait, and Kamchatka Strait. 
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Table 1. 26-year-mean fluxes across Aleutian Island cross-sections. Net volume flux 
(Net Vol.) is given in Sv, with the standard deviation shown below each 
value of net volume flux. Positive (northward) volume flux (Pos. Vol.) and 
negative (southward) volume flux (Neg. Vol.) values are also shown. The 
heat and freshwater flux values are shown in the same way. Heat flux 

3 

values are given in TW, and freshwater flux values are given in mSv (10 
m 3 /s). The locations of the sections are shown in Figure 3. 


Pass/Strait 

Net Vol. 
/s.d. 

Pos. Vol. 
/Neg. Vol. 

Net Heat 
/s.d. 

Pos. Heat 
/Neg. Heat 

Net FW 
/s.d. 

Pos. FW 
/Neg. FW 

Unimak 

0.069 

0.078 

0.096 

-0.026 

1.933 

2.406 

2.968 

-1.035 

4.110 

4.570 

5.492 

-1.381 

Akutan 

-0.001 

0.025 

0.011 

-0.012 

-0.209 

0.890 

0.305 

-0.514 

-0.013 

1.368 

0.589 

-0.601 

Umnak 

0.003 

0.024 

0.015 

-0.012 

-0.026 

0.856 

0.485 

-0.510 

0.195 

1.296 

0.797 

-0.603 

Strait H 

0.002 

0.003 

0.003 

-0.002 

0.058 

0.112 

0.124 

-0.065 

0.083 

0.156 

0.161 

-0.078 

Samalga 

0.041 

0.090 

0.344 

-0.303 

1.085 

3.013 

10.760 

-9.674 

2.355 

4.653 

15.145 

-12.790 

Strait G 

0.040 

0.100 

0.116 

-0.077 

1.131 

3.156 

3.721 

-2.590 

2.111 

5.256 

6.067 

-3.956 

Amukta 

1.601 

0.842 

1.753 

-0.152 

50.378 

25.783 

55.408 

-5.030 

64.783 

31.020 

71.284 

-6.501 

Seguam 

0.026 

0.059 

0.102 

-0.076 

0.724 

2.171 

3.441 

-2.717 

1.375 

2.975 

4.767 

-3.392 

Strait F 

0.006 

0.024 

0.029 

-0.023 

0.121 

0.854 

0.983 

-0.862 

0.306 

1.100 

1.315 

-1.009 

Strait E 

-0.004 

0.071 

0.058 

-0.062 

-0.426 

2.516 

1.944 

-2.370 

-0.002 

3.214 

2.659 

-2.661 

Strait D 

-0.006 

0.054 

0.046 

-0.053 

-0.367 

1.917 

1.579 

-1.946 

-0.196 

2.471 

2.117 

-2.313 

Amchitka 

1.890 

2.299 

2.962 

-1.072 

58.580 

64.668 

87.697 

-29.117 

50.394 

32.053 

63.806 

-13.412 

Strait C 

-0.136 

0.757 

0.464 

-0.600 

-4.194 

23.156 

14.344 

-18.537 

-2.881 

18.225 

12.410 

-15.291 

Buldir Pass 

1.216 

1.566 

2.178 

-0.962 

36.673 

47.016 

66.280 

-29.607 

22.635 

29.406 

45.088 

-22.453 

Strait B 

-0.003 

0.019 

0.009 

-0.012 

-0.081 

0.604 

0.302 

-0.383 

-0.092 

0.732 

0.332 

-0.424 

Near Strait 
1 

2.520 

4.418 

3.806 

-1.286 

70.282 

123.340 

106.062 

-35.780 

22.686 

49.537 

39.576 

-16.890 
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Pass/Strait 

Net Vol. 
/s.d. 

Pos. Vol. 
/Neg. Vol. 

Net Heat 
/s.d. 

Pos. Heat 
/Neg. Heat 

Net FW 
/s.d. 

Pos. FW 
/Neg. FW 

Near Strait 
2 

2.537 

6.607 

9.283 

-6.746 

66.455 

170.506 

238.023 

-171.568 

15.180 

54.676 

70.384 

-55.205 

Near Strait 
(total) 

5.057 

5.384 


136.737 

137.200 

344.085 

-207.348 

37.866 

48.680 

109.961 

-72.095 

Kamchatka 

-8.903 

4.849 

2.983 

-11.886 

-195.267 

107.142 

59.597 

-254.864 

-116.494 

37.752 

13.729 

-130.223 


1. Comparison of Observed and Modeled Volume Transports 

Multi-year continuous observations of the flow through the western passes are not 
available. Available observations include satellite-tracked drifters and short-term CTD 
measurements. Observations of the inflow through Buldir Pass were ~1 Sv during 
September 1992 (Stabeno and Reed 1993), based on CTD measurements (see Table 2 for 
a comparison of modeled and observed volume transport values). The model results 
show a long-tenn mean of 1.2 Sv [standard deviation (S.D.) = 1.566] net northward flow 
through Buldir Pass and show that September is the minimum of the annual cycle. 
Therefore, the observed value during September 1992 may be an underestimate of the 
long-term mean. 
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Table 2. Comparison of observed and modeled net volume transport values (Sv) for 
Buldir Pass, Near Strait, and Kamchatka Pass. Values for Buldir Pass and 
Near Strait are northward and values for Kamchatka Pass are southward. 
Model values are the mean of 26 years (1979-2004). Observed values are 
for time periods shown in parenthesis. 


Observation/Model 

Buldir Pass 

Near Strait 

Kamchatka Pass 

NAME Model 

1.2 (net), 2.2 
(northward), -1.0 
(southward) 

5.1 (net), 13.1 
(northward), -8.0 
(southward) 

-8.9 (net), 3.0 
(northward), -11.9 
(southward) 

Stabeno and Reed 
(1993) 

~1 

(Sept. 1992) 



Favorite (1974) 


10 

(1955-1974, 

discontinuous) 


Reed and Stabeno 
(1993) 


~5 

(Sept. 1992) 


Stabeno and Reed 
(1992) 


~ 3 

(Aug. 1991) 

-6.8 

(Aug. 1991) 

Verkhunov and 
Tkachenko (1992) 



~-ll 

(Oct.-Nov. 1990) 

-6 to -7 

(Apr-May 1990) 

Panteleev et al. 
(2006) 



-24 

(1932-2004, 

discontinuous) 


Early observations of the northward flow through Near Strait range from 3.7 to 26 
Sv (Arsen’ev 1967; Ohtani 1970; Hughes et al. 1974). Favorite (1974) concluded that a 
total of 10 Sv was a typical value. However, the methods behind all of these observations 
cannot be considered precise. More recently, Reed and Stabeno (1993) collected 
observations of volume flux through Near Strait during September 1992. Based on CTD 
casts, they found ~5 Sv of Alaskan Stream inflow through Near Strait. This compares 
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well with the 26-year mean modeled net inflow of 5.057 Sv (a combination of the flow 
through sections labeled Near Strait 1 and Near Strait 2; S.D. = 5.4; Table 1; Figure 3). 
Based on satellite-tracked drifters released throughout the region, Stabeno and Reed 
(1992) noted an anomalous lack of inflow through Near Strait that began in summer 1990 
and persisted at least through fall 1991. Modeled time series of volume flux through 
Near Strait over 26 years (Figure 5) show seven time periods when the transport was 
near-zero or even reversed. These anomalies last from three months to almost two years. 
While the model results do not show a lack of inflow at exactly the same time period as 
Stabeno and Reed (1992), the processes underlying the anomaly appear to be the same. 
These processes will be discussed in the next section. Stabeno and Reed (1992) note that 
a strong and prolonged reduction of inflow into the Bering Sea would likely lead to a 
large cooling of subsurface waters. A reduction of inflow would also promote 
salinization of surface waters, thereby increasing the depth of wintertime convection. 
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Figure 5. Monthly mean time series of net volume transport (Sv; red line) through 
selected sections, as indicated. The black line represents a 13-month 
running mean and the green line indicates the overall mean. Means for each 
time series are shown to the right of each panel. Local maxima (red arrows) 
and minima/flow reversals (blue arrows) are indicated on the Near Strait 

time series. 


Kamchatka Strait is the westernmost pass and is the primary location of 
southward flow out of the Bering Sea. Estimates of the Sverdrup transports within the 
deep Bering Sea basin imply that the local mean wind stress curl accounts for a 
significant portion (about 6 Sv) of the Kamchatka Current (Bond et al. 1994), which 
flows southward out of the Bering Sea via Kamchatka Strait. Outflow through 
Kamchatka Strait was estimated at ~11 Sv based on observations during April-May 1990 
and at 6-7 Sv during October-November 1990 (Verkhunov and Tkachenko, 1992). 
Stabeno and Reed (1992) made observations in Kamchatka Strait that yielded an estimate 
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of 6.8 Sv flowing southward during August 1991. However, they believe that this value 
is not typical of the Kamchatka flow, which instead is normally higher and closer to —11 
Sv. Panteleev et al. (2006) reconstructed the Bering Sea circulation as a variational 
inverse of the hydrographic and atmospheric climatologies, transport estimates through 
the Bering Strait, and surface drifter data. Based on this work, they estimate that the 
summer transport through Kamchatka Strait is 24 Sv southward, which is approximately 
twice as large as previous estimates. Panteleev et al. (2006) speculate that the difference 
between their estimate and previous observations is due to underestimation of the 
barotropic velocity component in the traditional transport estimates. The model results 
during 1979-2004 give a value of 8.9 Sv (S.D. = 4.8) net southward flow, which is within 
the range of the observations (Table 2), though much lower than the estimate by 
Panteleev et al. (2006). A more detailed comparison would be required to resolve this 
discrepancy, including the effect of only summertime observational results throughout 
the Bering Sea in the estimate by Panteleev et al. (2006) and their representation of the 
flow through other passes. 

It is important to note that flow through these passes and straits is highly variable 

on time scales from weeks to years. Stabeno and Reed (1992) found a dramatic decrease 

in the inflow of Alaskan Stream water into the Bering Sea from summer 1990-fall 1991. 

At the time of their publication, this was completely unexpected, since all previously 

available measurements showed relatively constant values of inflow and outflow. A high 

degree of variation occurs in the modeled time series of volume flux through the western 

passes and straits. Several strong flow reversals occur over the course of the 26-year time 

series in Near Strait and Buldir Pass. Five flow reversals occur in Kamchatka Strait, up to 

a few Sv northward, during the 26-year simulation. The monthly mean volume flux time 

series through Near Strait and Kamchatka Strait are significantly correlated (at the 99% 

significance level). The correlation coefficient is -0.77 for the monthly means values and 

-0.79 when the annual cycle is removed from both time series. The correlation 

essentially means that when the inflow into the Bering Sea via Near Strait is strong, the 

outflow via Kamchatka Strait also tends to be strong. Also, when the flow through one of 

these straits is weak, the other would be expected to be weak. This correlation from 
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model results is in qualitative agreement with observations by Stabeno and Reed (1992), 
which showed a weakened Near Strait inflow (~3 Sv) coincident with a weakened 
Kamchatka Strait outflow (6-7 Sv) in August 1991. 

2. Modeled Heat Fluxes 

Twenty-six-year mean modeled heat fluxes across the Aleutian Island Passes are 
shown in Table 1. The reference temperature used to calculate the heat flux is the 
freezing temperature, which is based on salinity in this calculation. The highest net 
northward heat fluxes into the Bering Sea in descending order are through Near Strait, 
Amchitka Pass, Amukta Pass and Buldir Pass (Figure 6). The net monthly mean heat 
flux through Near Strait ranges from 203 TW southward to 483 TW northward. The 
range of variability is ~5.0 times the mean net heat flux (137 TW) at Near Strait. The net 
monthly mean heat flux through Amchitka Pass ranges from 96 TW southward to 287 
TW northward. The range of variability is —6.5 times the mean net heat flux (59 TW) at 
Amchitka Pass. The net monthly mean heat flux through Amukta Pass ranges from 6 TW 
southward to 137 TW northward. The range of variability is ~ 2.9 times the mean net 
heat flux (50 TW) at Amukta Pass. The net monthly mean heat flux through Buldir Pass 
ranges from 74 TW southward to 142 TW northward. The range of variability is ~ 5.8 
times the mean net heat flux (37 TW) at Buldir Pass. Southward heat flux occurs 
primarily through Kamchatka Strait with a southward mean of 255 TW and a northward 
mean of 60 TW, which gives a net mean heat flux of 195 TW southward. When the net 
mean heat flux is summed for all of the cross-sections shown in Table 1, the total net 
mean heat flux across the Aleutian Island Passes is 87 TW northward. 
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Near Strait Heat Flux (Tref=Tfreeze; positive direction is North) 
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Figure 6. Monthly mean (red line), 13-month running mean (black line) and overall 
mean (green line) heat flux through selected passes as labeled. Fleat flux is 
referenced to the freezing temperature, which is based on salinity. 

3. Modeled Freshwater Fluxes 

Freshwater flux, using a reference salinity of 33.8, was calculated across the 
sections shown in Figure 3. Almost 70% of the northward freshwater flux from the North 
Pacific into the Bering Sea occurs between Unimak and Amchitka Passes, even though 
only 37% of the volume transport inflow occurs there. The highest mean freshwater 
fluxes entering the Bering Sea occur through Amukta Pass, Amchitka Pass, Near Strait 
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and Buldir Pass, in descending order (Figure 7). The net monthly mean freshwater flux 
through Amukta Pass ranges from 3 mSv southward to 180 mSv northward. The range of 
variability is ~2.8 times the mean net heat flux (65 mSv) at Amukta Pass. The net 
monthly mean freshwater flux through Amchitka Pass ranges from 54 mSv southward to 
135 mSv northward. The range of variability is ~3.8 times the mean net heat flux (50 
mSv) at Amchitka Pass. The net monthly mean freshwater flux through Near Strait 
ranges from 110 mSv southward to 186 mSv northward. The range of variability is ~7.8 
times the mean net heat flux (38 mSv) at Near Strait. The net monthly mean freshwater 
flux through Buldir Pass ranges from 48 mSv southward to 102 mSv northward. The 
range of variability is ~6.5 times the mean net heat flux (23 mSv) at Buldir Pass. When 
the net mean freshwater flux is summed for all of the cross-sections shown in Table 1, the 
total net mean freshwater flux across the Aleutian Island Passes is 67 mSv northward. 


25 



Freshwater Flux (mSv) freshwater ^ ux ( m Sv) Freshwater Flux (mSv) Freshwater Flux (mSv 


Amukta Pass Freshwater Flux (Sref=33.8; positive direction is North) 

I 

ii .nil! 



1£ 1 i I I I I L 


> 


~ 64.78 


1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 
Amchitka Pass Freshwater Flux (Sret=33.8; positive direction is North) 



50.39 


Near Strait Freshwater Flux (Sref=33.8; positive direction is North) 


100 

0 


JL 


✓°\ /\ 


1— 1 —I—r 

A 




1 A 


* V 


' I ' ' ‘ W ' 

jm. /W 


-100t _ i- i ■ i ■ i - L- -i Hir-.i - 

1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 


• 37.87 


Buldir Pass Freshwater Flux (Sref=33.8: positive direction is North) 


100 



— -a 22.63 


1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 


monthly mean 


13-mo. running mean 


mean 


Figure 7. Monthly mean (red line), 13-month running mean (black line) and overall 
mean (green line) freshwater flux through selected passes as labeled. 
Freshwater flux is referenced to 33.8 psu. 


D. MESOSCALE EDDIES IN THE ALASKAN STREAM AND THE DEEP 
BASIN 

1. Eddies in the Deep Bering Sea Basin 

The Bering Sea circulation over the deep basin is cyclonic in the mean, as shown 

in Figure 4. However, animations of monthly mean model results over the 26-year 
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simulation show frequent eddy activity throughout this area. For instance, at least 14 
mesoscale eddies are present in the Bering Sea during June 1987, as shown in Figure 8 
(Clement Kinney and Maslowski 2008). Half of these are anticyclonic and the other half 
are cyclonic. Diameters of these eddies are 120 km and greater and velocities are up to 40 
cm/s. Lifetimes of these eddies are typically a few months. The 9-km horizontal 
resolution of the model makes it possible to resolve eddies with diameters as small as 36 
km, however, the smallest eddies are likely not resolved. Recalling that the Bering Sea 
has a Rossby radius of deformation of-12-20 km according to Chelton et al. (1998). 



Reference Vectors 
(an/s) 


Figure 8. June 1987 upper 100 m circulation (vectors) and total kinetic energy 

(shading). Magenta contour lines represent bathymetry (m). The red line 
indicates the position of a cross section shown in Figure 5. Every second 

vector is shown. 


Observations of eddies in the deep Bering Sea basin have been made by Cokelet 
et al. (1996) and Stabeno and Reed (1994), among others (e.g., Solomon and Ahlnas 
1978; Kinder et al. 1980; Paluszkiewicz and Niebauer 1984). Stabeno and Reed (1994) 
observed a large anticyclonic eddy west of Bowers Ridge that had a diameter of -200 km 
and velocities of 30-40 cm/s. This eddy is similar in size and velocity to those from the 
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model; diameters are 120 km and greater, with velocities up to 40 cm/s. Cokelet et al. 
(1996) also observed a strong anticyclonic eddy here and suggested that it may be a 
recurring feature. 

Schumacher and Stabeno (1994) suggest several possible mechanisms for the 
formation of eddies in this region, including instabilities, wind forcing, strong flows 
through the passes, and topographic interactions. Cokelet et al. (1996) suggest that 
instabilities along the Bering Slope and Kamchatka Current and interactions with 
canyons and embayments at the landward edge of these currents, as well as inflows 
through the Aleutian Island Passes, may be responsible for eddy generation. Stabeno and 
Reed (1994) observed several eddies and meanders within the Kamchatka Current and 
the Bering Sea Basin by utilizing satellite-tracked drifters. Several anticyclonic eddies 
were observed in the western side of the basin and resulted from the interaction of the 
Kamchatka Current with topographic features. Stabeno and Reed (1994) suspect that 
these features are semi-permanent, since they appeared in drifter trajectories from more 
than one year. Model results (Figure 8) show two anticyclonic eddies associated with the 
Kamchatka Current during June 1987. (A cyclonic eddy (part of a counter-rotating pair) 
is also present near the Kamchatka Current at that time.) 

A vertical section of potential temperature through the southernmost eddy (Figure 
9; cross-section location shown in Figure 8) shows the cold core of the eddy with 
temperatures less than 1°C from -25-250 m depth. The anticyclone reversed the flow of 
this boundary current to over 15 cm/s northward, with a (local) southerly flow -75 km 
further offshore due to the eddy. Temperatures were up to 2.5 °C colder than the 26-year 
mean within the core of the eddy. Typical June temperatures in this location are 1.5-3°C 
and the usual speed of the southward-flowing Kamchatka Current is up to 9 cm/s in the 
core (Figure 9). A shallower thermocline is also apparent when the eddy is present, as 
shown by the temperature difference Figure 9c (negative values just below the surface). 
The temperature signature of the modeled eddy is similar to observations by 
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Rogachev et al. (2007) of cold-core eddies within the Kamchatka Current. One of the 
observed eddies was located just south of Kamchatka Strait and displayed temperatures 
near 0°C at its core. 

The salinity section reveals a strong displacement of isohalines (> 200 m) across 
the eddy during June 1987, as compared to the 26-year mean (1979-2004; Figure 10a). 
The relatively fresh core of the anticyclone exhibited salinities of 32.5 as deep as 200 m. 
Typically, water with such a fresh signal is confined to the shelf (Figure 10b). Salinity 
differences (between June 1987 and the 26-year mean June) show that the core of the 
eddy was up to 0.75 psu fresher than average, as well as -0.15 psu fresher at the surface. 
Rogachev and Gorin (2004) observed salinities of < 33 psu in an eddy near Kamchatka 
Strait during April 1991. 

Animations (not shown) of monthly mean velocity and sea surface height 
anomaly show that over the 26-year simulation (1979-2004), the Kamchatka Current was 
continuously populated with meanders and southwestward-moving anticyclonic eddies, 
such as the one just described. The diameters of these eddies were 83-139 km and the 
lifetimes were typically a few months. This is corroborated by observations of eddies 
with similar trajectories (Rogachev and Gorin 2004; Rogachev and Cannack 2002). 
Okkonen (1993) suggested that the anticyclonic eddies found in the western Bering Sea, 
such as the ones found within the path of the Kamchatka Current, may be related to 
planetary topographic waves. 
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Figure 9. Vertical cross-sections of potential temperature (°C; color shading) and 
velocity (cm/s; contours) along the red line in Figure 4 during (a) June 
1987, (b) 26-year mean June, and (c) the difference between these time 
periods. Solid contour lines indicate northward velocity and dashed lines 

indicate southward velocity. 
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Figure 10. Vertical cross-sections of salinity (color shading) and velocity (cm/s; 

contours) along the red line in Figure 4 during (a) June 1987, (b) 26-year 
mean June, and (c) difference between these time periods. Solid contour 
lines indicate northward velocity and dashed lines indicate southward 

velocity. 
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Eddies in the Alaskan Stream 


Maslowski et al. (2008a) discussed anticyclonic eddies that propagate along the 
path of the Alaskan Stream with an average diameter of 168 km and phase speed of ~2.3 
km day' 1 . These are similar to observations by Okkonen (1992), as well as those of 
Crawford et al. (2000) who detected the presence of six anticyclonic eddies over a period 
of seven years with an average diameter of 160 km and a mean phase speed of ~2.5 km 
day In the model simulations, these eddies did not reduce the strength of the Alaskan 
Stream, however there was an off-shore (or southward) shift in the velocity core. The 
mesoscale eddies were shown to have a strong effect on the volume and property fluxes 
through the eastern and central Aleutian passes (Maslowski et al. 2008a). Large increases 
in the northward volume, heat, and salt fluxes through Amukta and Amchitka passes 
occurred during the presence of an eddy in the Alaskan Stream (Maslowski et al. 2008a). 
O kk onen (1996) estimated that an eddy observed to separate from the Alaskan Stream 
south of Amchitka Pass represented 21% of the mean annual transport into the Bering 
Sea. 

Model results show that the larger and deeper passes of the western Aleutian 
Islands are also affected by these mesoscale features. The numerous eddies and 
meanders of the Bering Sea and within the Alaskan Stream play a critical role in 
detennining the flow through the western Aleutian passes and straits. As mentioned 
above, there were seven flow reversals in Near Strait throughout the 26-year simulation. 
Eddy activity in the vicinity of the strait has a strong impact on the flow through the 
strait. Figure 11 shows the local circulation in the upper 400 m during each of the seven 
flow reversals. There does not appear to be a consistent circulation regime associated 
with a flow reversal, instead the reversals occur due to the relative positions of multiple 
eddies and meanders both north and south of Near Strait. For example, during October 
1982 (Figure 11), mesoscale activity modified the flow through Near Strait: on the 
eastern side, a large anticyclone modified the path of the Alaskan Stream inflow and on 
the western side, anticyclonic flow is found within the strait, creating a strong southward 
flow. A total of 18.6 Sv flowed southward at this time, with 11 Sv flowing northward, 
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which equaled a net southward flux of 7.6 Sv. In March 1985, a large anticyclone has 
pushed part of the Alaskan Stream into Near Strait, causing a southward flow there. At 
the same time, the anticyclone is blocking northward flow through the strait. During May 
1989, an anticyclone is actually responsible for a northward flow (10.5 Sv) in the eastern 
part of Near Strait, however a stronger southward flow (14.3 Sv) exists across the western 
part of the strait (likely due to a cyclone to the north), which gives a net southward 
transport of 3.8 Sv. Anticyclonic blocking occurs again in July 1992, March 1995, and 
May 2002. August 1998 shows a strong southward flow in the central part of Near Strait, 
which originates in the Bering Sea, due to local cyclonicity. 
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Figure 11. Monthly mean total kinetic energy (shading; cm /s ) and velocity vectors 
(cm/s) in the upper 400 m during each of the flow reversals as shown in 
Figure 5. The year and month are indicated above each panel. The red line 
indicates the position of Near Strait and the green contour lines indicate 

bathymetry. 
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Four local maxima are indentified in the time series of volume transport through 
Near Strait (Figure 5). Similar to the flow reversals, these maxima are strongly affected 
by eddies and meanders both north and south of the strait. During June 1983, the 
northern part of a large, irregularly shaped anticyclone reached into the strait, thereby 
creating a strong net northward flow (17.7 Sv; Figure 12). Note that this same 
anticyclone contributed to a flow reversal 8 months earlier by blocking northward flow in 
the eastern part of the strait (Figure 11). The October 1988 circulation pattern is quite 
complex, however a strong northward flow in the central part of the strait appears to be 
related to an anticyclonic meander of the Alaskan Stream. January 1994 and March 1996 


both show upstream anticyclones that are pushing water northward through the eastern 


part of Near Strait. 



— 1.00 
—» 10.0 
-a- 20.0 

2 2 

Figure 12. Monthly mean total kinetic energy (shading; cm /s ) and velocity vectors 
(cm/s) in the upper 400 m during each of the flow maxima as shown in 
Figure 5. The year and month are indicated above each panel. The red line 
indicates the position of Near Strait and the green contour lines indicate 

bathymetry. 
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E. CONCLUSIONS 


High-resolution modeling of the Bering Sea basin has been limited. Hermann et 
al. (2002) presented results from their regional, eddy-resolving model for 1995 and 1997. 
However, their limited domain (the southeast Bering Sea) did not include the Bering Sea 
basin. Instead, they focused on comparing hindcast results from 1995 and 1997, 
primarily of the eastern portion of the Bering Slope Current and the ANSC. Model 
results were also compared with Eulerian and Lagrangian observations. 

Previous modeling work by Overland et al. (1994) focused on a limited domain: 
from 166°W to 157°E and from 46°N to the Bering shelf break. The primitive equation, 
three-layer hydrodynamic model had 1/8° resolution and excluded shallow regions <500 
m deep. The Alaskan Stream was specified as a boundary condition at inflow and 
outflow points with a constant volume transport of 15 Sv. This approach allowed 
Overland et al. (1994) to realize general circulations features (e.g., complex cyclonic flow 
in the Bering Sea basin) similar to observations. However, they were not able to simulate 
the meanders in the Kamchatka Current that were identified by Stabeno and Reed (1994), 
nor were they able to show enough interannual variability in the flow through Near Strait 
to account for observations by Stabeno and Reed (1992) or Reed et al. (1993). The 
prescribed boundary condition for the Alaskan Stream and the climatological atmospheric 
forcing that was used may have prevented Overland et al. (1994) from simulating realistic 
variability in the flow through Near Strait. Although early observations (e.g., Reed 1984; 
Stabeno and Reed 1991) showed little eddy energy of the Alaskan Stream where it 
entered the Overland et al. (1994) domain (near 166°W), more recent observations are in 
stark contrast. For example, Crawford et al. (2000) observed 6 multi-year anticyclonic 
meanders and eddies in the Alaskan Stream using data from the TOPEX/Poseidon 
altimeter between 1992 and 1998. All of these mesoscale features formed east of 166°W, 
with some forming east of Kodiak Island. The range of variability in the Near Strait 
inflow from Overland et al. (1994) was between ~5-7 Sv (see Figure 7 in Overland et al. 
1994) over a period of ten years (sampled daily), whereas NAME simulates variability 
between -9.5 and 18.3 Sv over a period of 26 years (sampled monthly; Figure 5). 
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Results presented herein are the first from a high-resolution, multi-decadal 
simulation of the entire Bering Sea. Utilizing a large pan-Arctic domain with a non- 
prescribed Alaskan Stream and Aleutian throughflow allows reproduction of the observed 
meanders and eddies in the Alaskan Stream and Kamchatka Current (Figure 8). Strong 
interannual variability in flow through Near Strait (Figure 5) was simulated by the model, 
which captures (and exceeds) the range of variability that has been observed (Stabeno 
and Reed 1992; Reed et al. 1993). The modeled net northward flow through Near Strait 
approached zero and reversed 7 times from 1979-2006 for time periods of three months 
to two years. Therefore, it is believed that short-tenn observations (months to years) may 
not be representative of the actual mean flow. This indicates a need for continuous 
monitoring of the flow through Near Strait. 

During the 26-year simulation, meanders and eddies are continuously present in 
the Kamchatka Current, as well as elsewhere throughout the Bering Sea. This is 
consistent with results from Cokelet and Stabeno (1997), which show that the 
background flow in the interior of the Bering Sea is dwarfed by the energetic eddies 
which populate the region. 
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V. SHELF-BASIN EXCHANGE AND OUTER SHELF DYNAMICS 


A. INTRODUCTION 

The research presented in this chapter seeks to identify and discuss the regions 
and processes of shelf-basin exchange within the Bering Sea. It is also an attempt to 
provide a quantitative measure of mean volume and property fluxes across the outer 
shelf, defined as a region between 50-m and 200-m isobaths. This exchange is important 
to biological productivity along the slope and on the Bering Sea shelf because the deep 
basin is the main source of nutrients for these locations. In addition, the downstream 
biological production within the Chukchi and Beaufort seas is strongly dependent upon 
the influx of nutrients from the deep Bering Sea via Bering Strait (Walsh et al. 1989; 
Codispoti et al. 2005; Walsh et al. 2004; Walsh et al. 2005). Without this nutrient source, 
biological production on the Bering Sea shelf and, consequently, in the Western Arctic 
would be much lower. Stabeno et al. (1999) state a need to identify and understand the 
mechanisms, which provide nutrients to the euphotic zone of the Bering Sea shelf. 
Model results presented herein will describe some of these mechanisms. 

A description of the general circulation from the long-term-mean point of view 
based on model results is presented in Section B. Sections C-E present analyses of 
volume, heat, salt, and freshwater transports through a number of cross-sections spanning 
the entire outer shelf, which is an attempt to determine the important pathways of shelf- 
basin exchange and describe the outer shelf dynamics. Section F contains a description 
of the mean state and variability of the Bering Slope Current and also includes a 
discussion of the eddy kinetic energy. Section G contains analysis of the exchanges 
through Zhemchug Canyon and the effects of eddy interactions along the slope and 
within canyons on shelf-basin exchange. Finally, the results of the study are discussed 
and major findings are indentified in Section H. This chapter was published by Clement 
Kinney et al. (2009). 
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GENERAL CIRCULATION 


Twenty-six-year (1979-2004) mean modeled circulation and TKE (0-220 m) are 
shown in Figure 13. The Alaskan Stream to the south of the Aleutian Island Chain and 
Anadyr Current (between Siberia and St. Lawrence Island) reach mean speeds in excess 
of 35 cm/s. The Bering Slope Current (BSC) is characterized by speeds up to ~12 cm/s 
in the long-tenn mean. The Bering Sea middle shelf domain ranges between < 1-7 cm/s, 
with higher speeds located north of the Zhemchug Canyon region (~170-175°W). 



—> 10 cm/s 

Figure 13. 26-year (1979-2004) mean 0-220m circulation and total kinetic energy. 

Every second vector is shown. 
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Previous modeling and observational studies (Kinder et al. 1986, see Figure 5; 
Schuert and Walsh 1993) mention that the BSC splits into two branches as it nears the 
coast of Siberia near Cape Navarin. One branch was believed to turn north and circulate 
around of the Gulf of Anadyr, while the other was thought to turn southward and 
eventually, form part of the Kamchatka Current. However, little direct observational 
evidence is available to substantiate this idea. Model results show that the vast majority 
of the BSC follows the bathymetry and turns south as it nears Cape Navarin (Figure 13; 
see Figure 14 for geographic place names). Very little to none of the water from the BSC 
turns northward upon nearing the Siberian coast. There is also evidence of a narrow 
coastal current flowing northward around the cape; however this northward flowing 
water is not derived from BSC water. Based on their observations, Stabeno et al. (1999) 
propose a schematic circulation in agreement with model results presented here, showing 
the majority of BSC water turning southward near Cape Navarin. 

In an effort to visually simplify the mean circulation, a schematic circulation 
based on Figure 13 was constructed (Figure 14). This is an update to a previously 
published schematic (Figure 1 in Clement et al. 2005) and the major difference is the 
above-mentioned southward turn of the Bering Slope Current as opposed to a north-south 
split. The current schematic shows the generally northwestward flow across the shelf, as 
well as the Anadyr Current and Alaska Coastal Current (ACC). The Alaskan Stream is 
strongest with significant throughflow (primarily northward into the Bering Sea) 
occurring through Unimak, Amukta and Amchitka passes. 
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Figure 14. Schematic circulation (0-220m) and total kinetic energy based on Figure 
13. The locations of model cross-sections and names are included. 
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C. VOLUME TRANSPORT 

We have created a series of cross-sections, which approximate the 50- and 200-m 

isobaths and several sections connecting the two isobaths. Monthly mean values of 

volume, heat, salt and freshwater flux across each section were calculated for the years 

1979-2004. Table 3 shows the 26-year mean values of these fluxes for each cross- 

section. The section locations are shown and labeled in Figure 14. Figure 15a shows the 

26-year mean value of volume transport (Sv) across each section with an arrow indicating 

the mean direction of flow. Arrows are scaled as a percentage of the largest cross- 

sectional value shown in the figure (e.g., a volume transport of 0.5 Sv would have an 
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arrow twice as wide as a section with volume transport of 0.25 Sv). The net volume 
transport is positive (on-shelf) in the central and eastern part of the Bering Sea along the 
200-m isobath. It is highest through the section north of Unimak Pass (AS 1). In contrast, 
the westernmost sections (AS7-AS9) have relatively small negative (off-shelf) mean 
volume transport values. This is largely due to the separation of the BSC from the shelf- 
break west of AS6, which removes the source of large (>0.05 Sv) on-shelf transports. 
The cross-shelf sections (CS1-CS6) increase in volume transport from east to west, as 
upwelled water is accumulated and integrated into a generally northwestward flow along 
the outer shelf. The sharpest increase between consecutive cross-shelf sections occurs 
between CS3 and CS4 (increase of 0.224 Sv) over a distance of approximately 230 km. 
This increase in shelf transport is due to water being moved on-shelf through the 
Zhemchug Canyon (AS4-AS6) to the south. Sections along the 50-m isobath have very 
small values of net volume transport, except for section AN6, which has a value of 0.528 
Sv or -83% of the total net northward volume transport across the 50-m isobath. 
Bathymetric steering is the main cause of the low volume transports through sections 
AN1-AN5, since the 50-m isobath directs water toward AN6 and eventually through 
Anadyr Strait, which is a deeper pass than the eastern Shpanberg Strait. 
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Table 3. 26-yr mean volume transport (Sv), salt flux (million kg/s), freshwater flux 

(mSv; relative to 34.8), and heat transport (TW; relative to -0.1°C). 
Numbers in parenthesis represent standard deviation. 
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Figure 15. 26-year (1979-2004) volume transport (Sv; a) and heat flux (TW; relative 
to -0.1 °C; b) across various sections. Arrows indicate net direction (positive 
is North or West) and are scaled relative to the largest value in each figure. 
The shading indicates depth (m). The 50 m and 200 m isobaths are shown 
as dotted lines. The mean flux divergence for each polygon is shown in 

italics and is circled. 
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D. HEAT TRANSPORT 


The heat flux referenced to -0.1 °C was calculated through the sections shown in 
Figure 14. Pathways of heat transport onto and across the Bering Sea shelf are somewhat 
different from those of volume transport (Figure 15b). This is due to differing pathways 
of warmer (less dense) water from colder (more dense) water. The majority of heat flux 
onto the shelf occurs through section AS 1 and is primarily water that has recently entered 
the Bering Sea via eastern Aleutian Island Passes and is relatively wann. Differences 
between volume and heat flux are also due to the fact that exchanges between the ocean 
and atmosphere change the amount of heat in the ocean over time, while the volume 
remains relatively constant. The largest on-shelf oceanic heat flux is through Zhemchug 
Canyon (a total of 6.1 TW through sections AS4-AS6), with the second largest on-shelf 
heat flux through Bering Canyon (5.1 TW, section AS1). In Figure 15b and Figure 16ab, 
several polygons are created by the cross sections. The mean values of heat transport, 
salt flux and freshwater flux across each section for each polygon are summed, such that 
means directed into/outside the polygon are positive/negative. The resulting sum (the 
flux divergence) is shown circled and in italics inside each respective polygon. The flux 
divergence values shown in Figure 15b indicate that a large amount of heat (2.5 TW) is 
lost to the atmosphere in the polygon northwest of Zhemchug Canyon that is bounded by 
CS3, AS4, AS5, CS4, and AN4. Lower heat flux divergence values occur in the 
easternmost and westernmost polygons, and also in the polygon north of the 50-m 
isobath. 
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Figure 16. 26-year (1979-2004) mean salt flux (million kg/s; a) and freshwater flux 
(mSv; relative to S=34.8; b) across various sections. Arrows indicate net 
direction (positive is North or West) and are scaled relative to the largest 
value in each figure. The shading indicates depth (m). The 50 m and 200 m 
isobaths are shown as dotted lines. The mean flux divergence for each 
polygon is shown in italics and is circled. 


A net loss of heat from the ocean to the atmosphere occurs as water generally 
moves northward throughout the Bering Sea toward the Chukchi Sea and Arctic Ocean. 
A total of 16.015 TW is advected onto the shelf across the 200-m isobath (Table 4). This 
calculation was made by summing the 26-year mean heat flux across the sections that lie 
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along the 200-m isobath (AS1-AS9; shown in Figures 14 and 15b). Approximately 75% 
of this heat is lost to the atmosphere between the 200- and 50-m isobaths. (The sections 
used in calculating the total heat flux across the 50-m isobath were AN1-AN6). A 
smaller loss occurs between the 50-m isobath and Anadyr and Shpanberg straits. 
Eventually, the net oceanic heat flux into the Chukchi Sea through Bering Strait is ~ 2.4 
TW or 14.9% of the 200-m isobath value. This implies that an increase in Bering Sea heat 
content could have a large impact on the Arctic Ocean, especially in light of the recent 
decline in sea ice. 


Table 4. Net heat transport (TW) across various locations in the Bering Sea. 


Location 

Heat transport (TW) 

Percentage of 200-m isobath 
value 

200-m isobath 

16.015 

100% 

50-m isobath 

3.934 

24.60% 

Anadyr + Shpanberg straits 

2.857 

17.80% 

Bering Strait 

2.383 

14.90% 


E. SALT AND FRESHWATER TRANSPORT 

The mean salt flux through the Bering Sea sections (Figure 16a) closely resembles 
the mean volume transport in a relative sense, due to the strong dependence of salt flux 
on volume flux. Significant on-shelf salt fluxes occur between Bering Canyon and 
Zhemchug Canyon and, as a result, cross-shelf fluxes (CS1-CS5) on the outer shelf 
increase as the water moves northwestward across the shelf. The vast majority of salt 
transport on the outer shelf (17 million kg/s or 83% of the total) eventually crosses the 
50-m isobath at AN6 near Anadyr Strait, with very little northward salt transport 
northward across the central and eastern sections. This is primarily due to bathymetric 
steering, as mentioned above in Section C; the 50-m isobath directs the shelf flow toward 
AN6 and eventually northward through Anadyr Strait. 
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The mean freshwater flux (referenced to S=34.8) through the Bering Sea sections 
shows that the primary freshwater route across the 200-m isobath is through the 
easternmost section (AS1; Figure 16b). This section accounts for 31% of the total net 
northward freshwater transport along the 200-m isobath. Similar to volume and salt 
transport, the total freshwater flux accumulates as water moves northwestward across the 
shelf through sections CS1-CS5. Then, the majority of freshwater turns northeastward 
and flows through Anadyr Strait. Figure 16ab shows that the polygon in the Gulf of 
Andyr (bounded by AN6, CS5, AS8, CS6, and AN7) has the highest values for salt and 
freshwater flux divergence. This is likely related to the large amount of sea ice produced 
and melted in this region, along with sea ice export/import into the polygon from 
surrounding areas. 

F. BERING SLOPE CURRENT 

The Bering Slope Current extends from the eastern Aleutian Islands along the 
shelf-break toward Cape Navarin. It is a relatively broad, northwestward flow and, more 
recently, has been described as a system of eddies (Okkonen 1993). Eddies ranging in 
size from 90 to 325 km are found to populate the Bering Sea basin, especially along the 
downstream portion of the Bering Slope Current to the northwest (see previous chapter). 
Several cross-sections across the Bering Slope Current and across several canyons that lie 
along the Bering Sea shelf-break were created to examine the flow and property fluxes in 
this region (locations of these cross-sections are shown in Figure 17). Figure 18 shows 
time series of monthly mean volume transport, heat flux, and section mean temperature 
through sections BSC1 and BSC4 from model results. The long-term mean volume 
transport for the BSC is 2-3.5 Sv. Flow reversals occur periodically at BSC1, however at 
the westermnost section (BSC4) reversals occur multiple times per year, likely due to 
passing eddies. In fact, monthly mean volume transport at BSC4 ranges between -15 to 
26 Sv. 
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Figure 17. Locations of Bering Slope Current and canyon cross-sections (solid black 
lines). The location of the 200-m isobath is shown as a thick dashed line. 
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Figure 18. Time series of monthly mean volume transport (a,d), heat flux (b,e) and 
mean temperature (c,f) across BSC1 (top 3 panels) and BSC4 (lower 3 
panels). Heat flux is referenced to -0.1°C. The red line represents monthly 
mean values, while the thick black line represents a 13-month running mean 
and the horizontal line represents the 26-year mean. Section locations are 

shown in Figure 17. 


The BSC carries a significant amount of heat (42-64 TW; Figure 18b,e), of which 
a portion is transported onto the shelf across the slope, as shown in Figure 15b. Volume 
transport is strongly correlated with heat transport at BSC1 and BSC4 (r = 0.99 for 
monthly mean time series at each section); however, the low-frequency (smoothed) 
section mean temperature shows a different, somewhat decadal, signal. Section mean 
temperature was generally above average from 1979-1989 at BSC1 and BSC4. It then 
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dropped below the 26-yr mean from 1989 to 1998/1999 and has recently been above 
average. Linkages to large-scale weather indices are examined and discussed in Chapter 
VIII. 


Correlation coefficients for volume transport, heat flux, and mean temperature 
between the four BSC sections are shown in Table 5. A significant (at the 95% level), 
but weak correlation exists for both volume transport and heat flux between BSC1 and 
BSC2, and these correlations become even less as fluxes across BSC1 are correlated with 
fluxes across sections further west. In fact, the correlations between BSC1 and BSC4 for 
both volume transport and heat flux are not significant (at the 95% level). This is possibly 
due to the increased eddy activity in the west, which acts to reduce the correlation with 
the upstream portion of the current. Another possible explanation for the reduced 
correlation between time series at eastern and western sections is that sometimes the BSC 
separates from the shelf break. This separation has been seen in animations of sea surface 
height, and also velocity throughout 26 years of model results. Over time scales of a few 
months, the position of the BSC is highly variable, especially in the western Bering Sea, 
and appears to be related to the fonnation of eddies just south of the shelf break. This 
separation from the shelf break is suggested in Figures 13 and 14, where TKE is high (up 
to 100 cm 2 s' 2 ) near the junction of AS3-AS4 and near AS6-AS7. If the BSC separates 
from the shelf break, conditions at BSC4 are not likely to be similar to conditions at 
BSC1. Mean temperature is strongly correlated between all BSC sections; however, it is 
important to note that there is no lag time. This, along with the lack of a strong 
correlation in volume transport, indicates that the temperature correlation is due to larger- 
scale basin changes in temperature, instead of an advected signal. 
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Table 5. Correlation coefficients for volume transport, heat flux, and section mean 
temperature between Bering Slope Current sections. The number in 
parenthesis is the lag time in months. The first column (Mo. Mean) for 
each parameter is the correlation for the monthly mean values. The second 
column (Annual Cycle Removed) is the correlation after removing the 
annual cycle. All correlations are significant at the 95% level except for 

those marked with a *. 


Section 

1 

Section 

2 

Volume transport 

Heat Flux 

Mean 

Temperature 

Mo. 

Mean 

Annual 

Cycle 

Removed 

Mo. 

Mean 

Annual 

Cycle 

Removed 

Mo. 

Mean 

Annual 

Cycle 

Removed 

BSC1 

BSC2 

0.32 (0) 

0.31 (4) 

0.25 (0) 

0.31 (4) 

0.82 (0) 

0.66 (0) 

BSC1 

BSC3 

0.22 (0) 

0.16(0) 

0.17(0) 

0.13 (0) 

0.80 (0) 

0.58 (0) 

BSC1 

BSC4 

0.08* (0) 

0.14(7) 

0.08* (0) 

0.14(7) 

0.75 (0) 

0.45 (0) 

BSC2 

BSC3 

0.32 (0) 

0.22 (0) 

0.29 (0) 

0.20 (0) 

0.83 (0) 

0.68 (0) 

BSC2 

BSC4 

0.23 (1) 

0.20(1) 

0.25 (0) 

0.23 (1) 

0.75 (0) 

0.51 (0) 

BSC3 

BSC4 

0.23 (5) 

0.36 (5) 

0.27 (5) 

0.38 (5) 

0.84 (0) 

0.65 (0) 


A previous study (Clement et al. 2005) addressed the eddy kinetic energy (EKE) 
of the Bering Sea during 1987 utilizing model results from the same model described 
herein. Along the BSC, the EKE was 20-100 cm" s' southeast of Zhemchug Canyon, 
but was significantly higher (80-300 cm" s'") northwest of Zhemchug Canyon. This is in 
agreement with Okkonen (1993) who found increased eddy activity in the downstream 
(western) leg of the BSC as compared to the upstream leg. However, it is possible that 
EKE is even higher in the real ocean, since this model is not yet fully eddy-resolving at 9- 
km horizontal resolution (Maslowski et al. 2008b). 

G. EXCHANGE IN CANYONS AND THE IMPORTANCE OF EDDIES 

Figures 15 and 16 indicate that volume transport and property fluxes are locally 
enhanced near canyons along the Bering Sea slope: Bering Canyon, located on the north 
side of the Alaska Peninsula, Pribliof Canyon, located south of the Pribilof Islands, and 
Zhemchug Canyon, located along the shelf break in the central Bering Sea between 
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171°Wand 176°W (see Figure 17 for canyon locations). Circulation and property 
exchange associated with these canyons is examined in greater detail below. 

Twenty-six-year (1979-2004) mean values of volume, heat, salt and freshwater 

transport through each canyon section are shown in Table 6. In the mean sense, the 

transport of water and all properties is northward through all of the canyons. However, 

reversals do occur and may sometimes even show up at the monthly mean time-scale 

(figures not shown). Because of the relatively strong on-shelf transports of volume, heat, 

and salt through section C6 (hereafter Zhemchug Canyon), the relationship between 

conditions in the BSC and exchange through Zhemchug Canyon are examined. Figure 19 

shows heat and salt flux anomalies through sections BSC3 and C6 (perpendicular 

sections whose locations are shown in Figure 17). The two smoothed time series of heat 

flux anomaly are somewhat correlated, with C6 lagging BSC3 by four months (r=0.51, 

significant at the 90% level). The salt flux anomaly time series are slightly less 

correlated (r=0.44, significant at the 90% level), again with a four-month lag. The salt 

flux anomaly is highly variable and peaks at over 10 million kg/s during November 1993 

(indicated with an open circle on Figure 20d). Based on Figure 20, the increase in salt 

transported northward onto the shelf in November 1993 is associated with a cyclonic 

eddy located just south of the canyon section. The eddy has a diameter of 145 km and 

persists for a period of ~3 months. In addition, three other peaks in the heat and salt flux 

anomalies are associated with cyclonic eddies (Figure 20). These peaks are represented 

with stars on the time series in Figure 20cd. A vertical section along the Zhemchug 

Canyon (BSC3) during the maximum salt anomaly (Nov. 1993) shows the high salinity 

anomaly core sitting above the slope, with positive salinity anomaly extending onto the 

shelf (Figure 21a). A section across Zhemchug Canyon (Figure 21b) shows that the 

salinity anomaly is positive near the surface (up to 0.12 psu) and also near the bottom (up 

to 0.18 psu), with the flow directed on-shelf. These figures show how the eddy is 

responsible for upwelling relatively salty water onto the central Bering Sea shelf. 

Focusing on the deeper water (>100m; where salinity can be used as a proxy for 

nutrients), one can say that this eddy accounted for a total of 28.6 trillion kg of salt being 

upwelled onto the shelf through the Zhemchug Canyon section during its three-month 
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lifetime. If the entire water column is included, then 107 trillion kg of salt was 
transported onto the shelf through Zhemchug Canyon during the same three-month 
period. 


Table 6. Twenty-six-year mean (1979-2004) volume transport (Sv), salt flux 
(million kg/s), freshwater flux (mSv; relative to 34.8), and heat transport 
(TW; relative to -0.1°C) for canyon sections. Numbers in parenthesis 
represent standard deviation. 


Section 

Volume 

transport 

Salt 

flux 

Freshwater 

flux 

Heat 

flux 

Cl 

0.206 

(0.084) 

6.601 

(2.705) 

15.842 

(6.495) 

5.245 

(2.892) 

C2 

0.009 

(0.080) 

0.281 

(2.579) 

0.679 

(6.120) 

0.556 

(1.505) 

C3 

0.083 

(0.087) 

2.682 

(2.799) 

6.198 

(6.492) 

1.261 

(1.593) 

C4 

0.135 

(0.053) 

4.357 

(1.717) 

9.839 

(3.840) 

2.501 

(1.233) 

C5 

0.110 

(0.071) 

3.556 

(2.303) 

7.895 

(5.122) 

2.035 

(1.444) 

C6 

0.195 

(0.112) 

6.284 

(3.616) 

13.927 

(7.976) 

3.550 

(2.369) 

Cl 

0.055 

(0.056) 

1.782 

(1.812) 

3.960 

(3.988) 

1.004 

(1.010) 
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Figure 19. Time series of monthly mean heat (a,c) salt (b,d) flux anomaly across BSC3 
(a,b) and the Zhemchug Canyon section (C6; c,d) shown in Figure 17. The 
red line represents monthly mean values, while the black line represents a 
13-month running mean and the horizontal green line represents zero. The 
circles and squares represent the salt and heat flux maximum (November 
1993) and minimum (November 1995), respectively through Zhemchug 
Canyon. The stars represent local peaks in the heat and salt flux anomaly 
through Zhemchug Canyon associated with eddies shown in Figure 20. 
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Figure 20. The bottom water salinity anomaly and sea surface height anomaly contours 
(blue is positive and red is negative; contour interval of 5 cm) during (a) 
March 1982, (b) July 1987, (c) November 1993, and (d) May 2002 in the 
vicinity of Zhemchug Canyon. The vectors represent monthly mean 
velocity over the entire water column during each respective month. 
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Figure 21. The vertical section of velocity (contours; cm/s) and salinity anomaly 

(shading) (a) along the Zhemchug canyon (section BSC3), and (b) cross the 
Zhemchug canyon (section C6) during November 1993. 


The contribution of a transient motion (e.g., an eddy) to the time-averaged flux 
across a particular section can be quantified according to: 

F = | {uc^dA (2) 

where F is a temporal mean flux (26-yr annual cycle; 1979-2004), J (...) dA is the area 
integral across the section, u is the normal component of velocity and c is the quantity 
being advected. Anomalies from the temporal mean are defined as: 
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u =u-u; 


( 3 ) 


c = c — c 

Taking Equations (2) and (3) together gives: 

F = J (wc)<fA + c jjA (4) 

The perturbations, u' and c', should be derived from instantaneous measurements 
of u and c; or in the case of model results, calculated based on high frequency temporal 
snapshot results. To make a fair calculation of the eddy transports, it would be best to 
use means (for u and c) that are daily, if not over a shorter time period. However, the 
model results are not available in three dimensions as daily means, due to storage 
limitations. Therefore, the monthly means must be used for calculating perturbation 
values. For the means ( u, c ), the 26-year annual cycle monthly means (e.g., mean Jan. 
over 26 years) are used. Transient motions contribute to the mean flux if u' and c' are 
correlated. The correlation coefficient between u' and S' across Zhemchug Canyon 
section was 0.32, while the correlation coefficient between u' and T was 0.20. Both 
correlations were significant at the 99% level (p=.01). The fractional contribution of 
eddy transports to the mean flux, F, was small for both heat and salt flux (< 8%). It is 
believed that the weak correlations and small eddy contributions to the mean flux are 
largely due to the fact that monthly means are used for calculation of perturbation values, 
which smooth out the short-lived perturbations (Marble 2001). This calculation does not 
fully quantify the eddy contribution to the mean flux; however, it is the best that can be 
made with the available model output. It is still apparent that eddies are responsible for 
several peaks in the heat and salt flux anomalies as shown in Figures 19-21. 

Hovmuller plots of near-bottom salinity and temperature anomaly (after removing 
the annual cycle) for 1979-2004 show significant interannual and spatial variability along 
the 200-m isobath (Figure 22). Zhemchug Canyon experiences frequent events, likely 
due to strong eddy activity, in comparison to Bering Canyon, which is more consistent. 
Overall, the temperature and salinity anomalies tend to be the same sign (either negative 
or positive) along the entire length of the 200-m isobath; however, this is not always true. 
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For instance during 1987, temperature and salinity anomalies were negative in the west 
(0-400 km), with mainly positive anomalies in the east. When looking over longer time 
scales, the temperatures were wanner, especially in the west during the 1980s, followed 
by a drop in the early 1990s. Recently, however, temperatures along the slope have 
increased. Salinity anomalies show some temporal similarities with temperature 
anomalies. Decreased salinity anomalies were present over the slope during the 1990s 
and increased in 2001. However, unlike temperature anomalies, the salinity was 
relatively low in the 1980s. There is also a suggestion of propagating features in these 
figures, possibly planetary waves or eddy propagation, which would be a topic for a later 
paper. 
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Figure 22. Hovmuller plot of temperature anomaly (°C; a) and salinity anomaly (b) 

along the Bering Sea 200-m isobath. 
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H. CONCLUSIONS 


The flow across the Bering Sea shelf is up to 3 cm/s in the mean sense and is 
highest near 60°N between St. Matthew and Nunivak islands. This area of strengthened 
flow is centered on and follows isobaths between 50 and 75 m. As expected, the 
strongest flow on the shelf is Anadyr Current to the west, which carries water toward 
Anadyr Strait at a mean speed of approximately 3-10 cm/s, depending upon the location. 
Several regions along the slope (especially certain canyons) that are important in shelf 
basin exchange also show up in the long-term mean circulation scheme (Figure 13). 

In the model results, the BSC appears to be more a system of eddies rather than a 
continuous current, which emphasizes the need for a fully eddy-resolving, basin-wide 
model to represent its complex dynamics. Okkonen (1993) identified a representative 
eddy period as ~4 months. This time period is similar to modeled high frequency signals 
in the volume transport and heat flux at BSC4 (Figure 18d, e), as well as the heat and salt 
flux anomalies at BSC3 (Figure 19). The lack of correlation between various sections 
along the BSC, along with snapshots of SSH indicating frequent, multiple eddies (not 
shown), provide additional evidence for this conclusion. 

There is little evidence for a north-south splitting of the BSC near Cape Navarin 
(Figure 13). Instead the majority of the BSC turns southwestward upon nearing the coast 
in the long-term mean model results. There is evidence of a coastal current flowing 
northeastward around Cape Navarin, however it is not derived of BSC water. According 
to the model results the main source of water for the Gulf of Anadyr comes from the 
outer shelf flow, roughly between the 100 and 200-m isobaths. 

Paluszkiewicz and Niebauer (1984) and Okkonen et al. (2004) have suggested 
that topographic planetary waves generate eddies seen along the Bering slope. The 
modeled eddy described herein had a diameter of 145 km and lifetime of ~90 days. This 
is similar in size and duration to eddies observed by Paluszkiewicz and Niebauer (1984; 
140 km / 84 days) and Okkonen (1993; 184 km / 72 days). Frequent eddies, such as the 
one described herein, may be important for transporting salt onto the Bering Sea shelf. 
However, if three-dimensional daily mean fields had been available to correlate velocity 
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and salinity along the Zhemchug Canyon section, instead of monthly means, this might 
have allowed a more quantitative statement to be made. On-shelf transport appears to be 
enhanced by the presence of canyons along the Bering Sea slope (Figures 15a, 16). 
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VI. CIRCULATION AND ENERGETICS OF THE NORTHERN 

BERING SEA 


A. INTRODUCTION 

This chapter aims to quantify the long-tenn mean and variability of the circulation 
at time scales ranging from months to decades across the data-limited northern Bering 
Sea and through Bering Strait using the NPS pan-Arctic coupled ice-ocean model. One of 
the specific objectives in this study was to improve understanding of the flow and 
transports through Anadyr, Shpanberg, and Bering Straits (locations shown in Figure 23). 
Model results on the EKE are also provided to indicate regions of active horizontal and 
vertical mixing, which is a biologically important process. These two topics are 
geographically related within the northern Bering Sea and are part of a publication by 
Clement et al. (2005). 

B. VOLUME TRANSPORT THROUGH ANADYR AND SHPANBERG 

STRAITS 

In order to determine the characteristics of water flowing through Bering Strait, it 
is important to consider upstream conditions. Model results show that water from Anadyr 
Strait contributes 80% of the mean flow (23-year mean = 0.52 Sv) into the Chukchi Sea. 
The wider, but shallower Shpanberg Strait contributes only 0.13 Sv or 20% of the mean 
flow through Bering Strait and is primarily made up of relatively sluggish Alaska Coastal 
Water and Bering Shelf Water. It is worth noting that not all of the water flowing through 
Anadyr Strait is associated with the Anadyr Current. Instead some water moving 
northward on the northern Bering shelf also can go west around St. Lawrence Island 
through Anadyr Strait to reach Bering Strait. Overland and Roach (1987) used a two- 
dimensional barotropic model to investigate the Bering Strait transport and circulation 
patterns of the northern Bering Sea. They found that in their 18-km grid in the absence of 
wind, Anadyr Strait contributed 72% of the northward transport that would eventually 
enter the Chukchi Sea via Bering Strait. Historic current meter measurements during 


65 



summer and fall from 1985-1988 indicate that the split between Anadyr/Shpanberg 
straits is typically 30/70%, respectively (Coachman 1993). 

Since the majority of Bering Strait throughflow is associated with the Anadyr 
Current, it is important to observe the characteristics of this water. However, most of 
Anadyr Strait lies in Russian territorial waters and international research efforts have 
been significantly limited in this area. Additional scientific efforts and field observations 
including long-term moorings in the Gulf of Anadyr and on the western side of Bering 
Strait, would be critical for determining the characteristics of the Pacific inflow into the 
Chirikov Basin north of St. Lawrence Island and eventually into the Arctic Ocean 
through the Bering Strait. Such data would allow further validation and constraining of 
ocean models for climate study. 

C. ANOMALOUS CIRCULATION EVENT 

Time series of the volume transport through five cross-sections (locations shown 
in Figure 23) are shown in Figure 24. Over 23 years, the sections show seasonal, as well 
as interannual, variability in flow across the northern Bering shelf. The net 23-year mean 
volume transport across Bering Strait is 0.65 Sv, with 0.72 Sv moving northward and a 
slight but annually regular component (0.07 Sv) moving southward. An annual minimum 
of 0.48 Sv was reached in 1994 and a maximum of 0.78 Sv in 1979, thus giving an 
interannual variability of approximately +/-0.15 Sv. A similar estimate of +/-0.2 Sv was 
made by Coachman (1993) based on observational studies during the Inner Shelf Transfer 
and Recycling in the Bering and Chukchi Seas (ISHTAR) program. Net monthly mean 
variation in transport is greater than 1 Sv, with a low of 0.11 Sv in December 1997 and a 
high of 1.28 Sv in August 1999. Examination of instantaneous transports across Bering 
Strait (calculated every model time-step of 480 s) reveals even stronger wind-driven 
variability, from 2.7 Sv northward to 2.5 Sv southward (data not shown). 
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Figure 23. Bathymetry (m) of the northern Bering Sea and locations of cross-sections 
(thick black lines): BS (Bering Strait), AS (Anadyr Strait), SS (Shpanberg 
Strait), AC (Anadyr Current), and SL (St. Lawrence Island). The direction 
of positive flow is northward across BS and SL, northeastward across AS 
and SS, and eastward across AC. 
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Figure 24. Monthly mean volume transport over a 23-year time series (1979-2001). 

Positive fluxes represent flow to the north or east according to the model 
grid (see Figure 23), while negative fluxes represent flow to the south or 
west. The smoothed net flux (thick green line) is a 13-month running mean. 
The 23-yr mean is shown at the end of each time-series. 


68 




































































































































































With respect to other sections, SS, AC, and SL all had numerous negative 
monthly mean minima representative of flow reversals (Figure 24). The most notable 
flow reversal (up to 0.70 Sv) in the 23-year time series occurred along the AC section 
during November and December 2000. Flow continued to be westward across the 
northwestern Bering Sea until March 2001. This episode was also evident in the AS, SS, 
and SL sections, with greatly increased northwestward flow across the northern Bering 
shelf. This event also coincided with observations of reduced sea ice conditions during 
winter of 2000-2001 (Figure 25). Observations mainly derived from satellite data during 
February 1999 and 2001 show that ice concentrations were very low during 2001 as 
compared to a more typical year, 1999 (Figure 25a,b; Clement et al. 2004). The model 
shows a similar ice distribution in the two years, except for slightly lower ice 
concentrations (80-95%) instead of solid >95% ice cover reported during February 1999 
(Figure 25c,d). These areas of slightly reduced ice cover concentrations in the model 
might be representative of the realistic presence of leads not observed by the U.S. 
National Ice Center. The strong flow reversal, which peaked in November and December 
2000, and northwestward movement of sea ice occurred relatively quickly over the 
course of a few months during winter 2000-2001. The cause of this event is most likely 
wind forcing, as shown in Figure 26. While previous studies indicate that prevailing 
winds are out of the north to slightly northeast in the northern Bering Sea during winter 
(Muench and Ahlnas 1976; Pease 1980; Overland 1981), a very different scenario took 
place during winter 2000-2001. Figure 26 depicts the ECMWF wind forcing fields used 
in the model averaged over the 2-month peak in flow reversal (November and December 
2000) as compared to the more typical year (November and December 1998). In 1998, 
wind was out of the north with speeds of 3-4 m/s in the vicinity of Bering Strait and up to 
4.5 m/s south of St. Lawrence Island. In contrast, during 2000 much stronger wind speeds 
of 5-6.5 m/s occurred near Bering Strait, with speeds up to 9.3 m/s in the Gulf of Anadyr. 
The average magnitude difference between 1998 and 2000 winter winds (Figure 26C) 
was 3.0 m/s and the wind direction was shifted an average of 43° to a more easterly 
pattern in 2000. Observational measurements from this time period corroborate an 
unusual wind and resulting ice cover pattern in winter 2000-2001 (Clement et al. 2004). 
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The dramatic response in sea ice and ocean conditions to this weather pattern suggests 
that wind is the dominant forcing mechanism in the northern Bering Sea. The wind 
forcing drove an Ekman transport on the shelf that resulted in an Anadyr Current flow 
reversal, as well as greatly increased northward transport across Shpanberg Strait. 



Figure 25. Sea ice concentration on (A) February 19, 1999 and (B) February 19, 2001 
as determined using data obtained from the U.S. National Ice Center 
(adapted from Clement et ah, 2004). Model monthly mean sea ice 
concentration during (C) February 1999 and (D) February 2001. 
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Figure 26. Model wind forcing fields averaged over November and December (A) 
1998 and (B) 2000. (C) The magnitude difference (2000-1998) and both 
wind fields (1998 in black; 2000 in white). 

D. EDDY KINETIC ENERGY 

With a horizontal grid cell spacing of approximately 9 km, this model is able to 
resolve eddies with diameters as small as 36 km. As such, the eddy kinetic energy (EKE) 
distribution in the Bering Sea has been examined, as calculated from daily fluctuations 
( u’,v ’) referred to the 23-year mean. EKE is an important model parameter, which, in the 
physical world, commonly represents regions of active mixing and associated increased 
biological productivity. Realistic representation of EKE in a model is often challenging 
but when achieved can provide useful infonnation about the dynamics and productivity in 
a region. Observational estimates of EKE on the Bering shelf are not readily available, 
hence the following discussion serves to establish a quantitative, yet due to the need for 
even higher model spatial resolution, possibly incomplete reference. Results described 
below use either the surface level (0-5 m) or a subsurface level (20-26 m) as indicated. 
The year 1987 was chosen because it characterized an average transport through Bering 
Strait and it did not appear strongly biased toward any particular climate regime in the 
North Pacific or Arctic Ocean. The mean velocity components for 1979-2001, were 
calculated and then subtracted from the 1987 daily velocity values to obtain u’ and v’: 

u ’ = u d - u m , v , =v d -v m (5) 
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where u d and v,i represent the daily velocity components and u m and v m represent the long- 
terni (1979-2001) mean velocity components. From this the EKE can be calculated as: 


EKE = 


( U ~ + V~) 

2 


( 6 ) 


It is important to note that high EKE values indicate not only the presence of 
eddies, but also indicate any fluctuation around the mean velocity field. 

Seasonal changes in EKE in the surface layer are visible across the Bering shelf, 
with highest values in autumn and winter (Figure 27). Bering and Anadyr straits maintain 
relatively high EKE throughout the year, but peak during autumn. The area directly south 
of St. Lawrence Island is also a region of increased EKE, while north of the island in the 
Chirikov Basin a distinct decrease is observed. Below the surface at 20-26 m depth, there 
is a sharp decrease in EKE across the middle shelf (Figure 28). Again, the EKE in the 
northern straits and along the BSC remains elevated. 
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Figure 27. Seasonally average EKE at the surface (0-5 m) calculated from daily 1987 
snapshots (against the 23-year mean; 1979-2001): (A) Winter (J-F-M) 
average, (B) Spring (A-M-J) average, (C) Summer (J-A-S) average, and (D) 

Autumn (O-N-D) average. 
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Figure 28. Seasonally average EKE at 20-26 m calculated from daily 1987 snapshots 
(against the 23-year mean; 1979-2001): (A) Winter (J-F-M) average, (B) 
Spring (A-M-J) average, (C) Summer (J-A-S) average, and (D) Autumn (O- 

N-D) average. 
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E. CONCLUSIONS 


Model results indicate that water from Anadyr Strait contributes 80% of the mean 
flow into the Chukchi Sea via Bering Strait. The remainder of the northward flow 
through Bering Strait (20%) comes from the eastern side of the northern Bering Sea shelf 
via Shpanberg Strait. More observations of the flow through Anadyr Strait are needed, 
since it provides the majority of the Bering Strait throughflow; however, this is 
complicated by a political boundary (the Russian EEZ), which has historically restricted 
US access and limited international collaboration. 

Time-series of volume transport show that the model is able to represent 
anomalous events, such as major flow reversals, which are corroborated by observations. 
Wind appears to be the dominant force driving water and ice movement across the 
northern Bering Sea at synoptic to interannual time scales. In addition to wind forcing, it 
is also important to consider buoyancy forcing along the Alaska coast, however this 
model is not able to fully resolve the narrow Alaska Coastal Current. 

EKE fields show that the northern Bering Sea maintains year-round high energy 
and significant horizontal and vertical mixing, especially in Bering and Anadyr straits. 
The highest energy levels occur during autumn and winter throughout the shelf, 
coincident with observations of increased storminess during those seasons (Overland 
1981; Overland and Pease 1982; Stabeno et al. 1999). Possible effects of these EKE 
fields on biological processes are discussed in Chapter IX. 
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VII. THE FLOW THROUGH BERING STRAIT 


A. INTRODUCTION 

The Pacific Arctic Region spans the sub-Arctic Bering Sea northward through the 
Chukchi and Beaufort seas and the Arctic Ocean. A narrow passageway, the Bering 
Strait, connects the wide and shallow shelves of the Bering and Chukchi seas. The 
Bering Strait is the only Pacific connection to the Arctic Ocean. The narrow (~85 km 
wide) and shallow (~50 m deep) strait provides low-salinity and high-nutrient Pacific 
Water to the Chukchi Sea and the Arctic Ocean. Many global and regional models face 
challenges with resolving oceanic exchanges across this narrow and shallow strait, 
mainly due to the requirement of high spatial resolution and the associated high 
computational cost to resolve it. In fact, many coarse-resolution models either have a 
closed Bering Strait or use a prescribed boundary condition. For example, Goosse et al. 
(1997) demonstrated a significant improvement in modeled ocean dynamics in a coarse 
resolution (3°x3°) model with an opened Bering Strait. They also found that opening 
Bering Strait produced a more realistically positioned sea ice edge in the Bering Sea, 
because wann water was allowed to advect further north onto the Bering-Chukchi shelf. 
Arctic freshwater budgets were also improved, with increased freshwater storage in the 
Greenland and Norwegian Seas. 

Scientific access across Bering Strait has been restricted due to the political 
boundary between the United States and Russia. The Russian-US Convention line, 
dividing the Exclusive Economic Zones (EEZs) of the two countries, lies between two 
islands near the center of the strait: Ratmanova Island (part of Russia, also called Big 
Diomede in the U.S.) and Little Diomede Island (part of the U.S.). While U.S. research 
has maintained moorings in the Bering Strait almost continuously since 1990, only for 
limited portions of that time has U.S. access been granted to the western side of the strait. 

The main goal of this work is to analyze output on the Bering Strait throughflow 
from several state-of-the-art regional and global Arctic-focused models and to compare 
model results against the available limited observational data. The volume and property 
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fluxes will be analyzed over a long time series (up to 26 years depending on available 
results from individual models). In addition to interannual changes, annual cycles in 
these parameters will also be examined. Data are from moored instruments placed near¬ 
bottom in three point locations in the vicinity of the strait (Figure 29a). Both 
observations and modeling have their own limitations in Bering Strait. Numerical models 
are limited by relatively coarse resolution in the strait, errors in forcing and omitted 
processes (e.g., tides), whereas observational results are limited by spatial coverage 
across the strait, and lack of upper layer measurements. 
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IBCAO BESTMAS 



Figure 29. Bathymetry (m) in the vicinity of the Bering Strait (a). Depth contours are 
every 10 m from the International Bathymetric Chart of the Arctic Ocean 
(Jakobsson et al. 2000). Model bathymetry (m) from (b) BESTMAS, (c) 
ECC02, (d) NAME, (e) ORCA, and (f) PIOMAS. The approximate 
locations of the moored observations are indicated with black circles. The 
cross-sections across Bering Strait are shown as black lines in each model 
bathymetry figure. Different orientations in figures b-f are due to different 
coordinate systems used in each model. 
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B. MODEL DESCRIPTIONS 


In this section, important aspects of five global and regional sea ice-ocean 
coupled models are described. Results from each of these models were provided to the 
author in order to investigate Bering Strait inflow (Table 7). 


Table 7. Basin infonnation on the five models used in this study. 


Model 

Global/Regional 

Atmospheric Forcing 

Resolution in Bering 
Strait 

BESTMAS 

regional 

NCEP/NCAR reanalysis 

~7 km 

ECC02 


Japanese 25-year reanalysis 

~18 km 

NAME 

regional 

ECMWF reanalysis 

~9 km 

ORCA 

global 

DRAKKAR Forcing Set (DFS 
3.1) reanalysis 

~13 km 

PIOMAS 

regional 

NCEP/NCAR reanalysis 

~22 km 


1. Bering Ecosystem STudy Ice-Ocean Modeling and Assimilation 
System (BESTMAS) 

BESTMAS (Zhang et al. 2010) is based on the coupled Parallel Ocean and sea Ice 
Model (POIM) of Zhang and Rothrock (2003). The sea ice model is the multicategory 
thickness and enthalpy distribution (TED) sea ice model (Zhang and Rothrock 2001; 
Hibler 1980). It employs a teardrop viscous-plastic rheology (Zhang and Rothrock 2005), 
a mechanical redistribution function for ice ridging (Thorndike et al. 1975; Hibler 1980), 
and a LSR (line successive relaxation) dynamics model to solve the ice momentum 
equation (Zhang and Hibler 1997). The TED ice model also includes a snow thickness 
distribution model following Flato and Hibler (1995). The ocean model is based on the 
Parallel Ocean Program (POP) developed at Los Alamos National Laboratory (Smith et 
al. 1992; Dukowicz and Smith 1994). Given that tidal energy accounts for 60-90% of the 
total horizontal kinetic energy over the southeastern shelf region of the Bering Sea 
(Kinder and Schumacher 1981), tidal forcing arising from the eight primary constituents 
(M2, S2, N2, K2, Kl, 01, PI, and Ql) (Gill 1982) is incorporated into the POP ocean 
model. The tidal forcing consists of a tide generating potential with corrections due to 

both the earth tide and self-attraction and loading following Marchuk and Kagan (1989). 
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The model domain of BESTMAS covers the northern hemisphere north of 39°N. The 
BESTMAS finite-difference grid is based on a generalized orthogonal curvilinear 
coordinate system with a horizontal dimension of 600x300 grid points. The “north pole” 
of the model grid is placed in Alaska. Thus, BESTMAS has its highest horizontal 
resolution along the Alaskan coast and in the Bering, Chukchi, and Beaufort seas, with an 
average of about 7 km for the whole Bering Sea and 10 km for the combined Chukchi 
and Beaufort seas. There are 20 grid cells across Bering Strait (Figure 29b), which allows 
a good connection between the Bering Sea and the Arctic Ocean. The TED sea ice model 
has eight categories each for ice thickness, ice enthalpy, and snow depth. The centers of 
the eight ice thickness categories are 0, 0.38, 1.30, 3.07, 5.97, 10.24, 16.02, and 23.41 m. 
The POP ocean model has 30 vertical levels of varying thicknesses to resolve surface 
layers and bottom topography. The first 13 levels are in the upper 100 m and the upper 
six levels are each 5 m thick. The model bathymetry is obtained by merging the IBCAO 
(International Bathymetric Chart of the Arctic Ocean) dataset and the ETOP05 (Earth 
Topography Five Minute Gridded Elevation Data Set) dataset (see Holland 2000). 
BESTMAS is forced by daily NCEP/NCAR reanalysis (Kalnay et al. 1996) surface 
forcing fields. Model forcing also includes river runoff of freshwater in the Bering and 
Arctic seas. For the Bering Sea, monthly climatological runoffs of the Anadyr, Yukon, 
and Kuskokwim rivers are used (Zhang et al. 2010). For the Arctic Ocean, monthly 
climatological runoffs of the Pechora, Ob, Yenisei, Olenek, Yana, Indigirka, Kolyma, 
Mackenzie, Dvina, Lena, Khatanga, Taimyra, and Piasina rivers are from the Alfred 
Wegener Institute (Prange and Lohmann 2004). Although BESTMAS has a large model 
domain that includes the Arctic and the North Pacific, realistic lateral open boundary 
conditions are still necessary to create the right water masses and fluxes. The POP ocean 
model has been further modified to incorporate open boundary conditions so that 
BESTMAS is able to be one-way nested to a lower resolution, global POIM (Zhang 
2005). Monthly mean open boundary conditions of ocean temperature, salinity, and sea 
surface height from the global POIM are imposed at the southern boundaries along 39°N. 
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2. Estimating the Circulation and Climate of the Ocean, Phase II 
(ECC02) 

The ECC02 regional Arctic Ocean solution uses a configuration of the 
Massachusetts Institute of Technology general circulation model (MITgcm; Marshall et 
al. 1997; Losch et al. 2010; Nguyen et al. 2010). The domain boundaries are at -55° 
North in both the Atlantic and Pacific sectors. These boundaries coincide with grid cells 
in a global, cubed-sphere configuration of the MITgcm (Menemenlis et al. 2005). 

The grid covering the Arctic domain is locally orthogonal with horizontal grid 
spacing of approximately 18 km. There are 50 vertical levels ranging in thickness from 
10 m near the surface to approximately 450 m at a maximum model depth of 6150 m. 
The model employs the rescaled vertical coordinate "z*" of Adcroft and Campin (2004) 
and the partial-cell formulation of Adcroft et al. (1997), which pennits accurate 
representation of the bathymetry. Bathymetry is from the S2004 (W. Smith, 2010, 
personal communication) blend of the Smith and Sandwell (1997) and the General 
Bathymetric Charts of the Oceans (GEBCO) one arc-minute bathymetric grid. The non¬ 
linear equation of state of Jackett and McDougall (1995) is used. Vertical mixing follows 
Large et al. (1994). A 7th-order mono tonicity-preserving advection scheme Daru and 
Tenaud (2004) is employed and there is no explicit horizontal diffusivity. Horizontal 
viscosity follows Leith (1996) but is modified to sense the divergent flow (Fox-Kemper 
and Menemenlis 2008). 

The ocean model is coupled to the MITgcm sea ice model described in Losch et 
al. (2010). Ice mechanics follow a viscous-plastic rheology and the ice momentum 
equations are solved numerically using the line-successive-over-relaxation solver of 
Zhang and Hibler (1997). Ice thennodynamics use a zero-heat-capacity formulation and 
seven thickness categories, equally distributed between zero to twice the mean ice 
thickness in each grid cell. Ice dynamics use a 2-category thickness with one for open 
water and one for ice. Salt rejected during ice formation is treated using a sub-grid-scale 
salt-plume parametrization described in Nguyen et al. (2009). The model includes 
prognostic variables for snow thickness and for sea ice salinity. 
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Initial and lateral boundary conditions come from the globally optimized ECC02 
solution (Menemenlis et al. 2008). Surface atmospheric forcing fields are from the 
Japanese 25-year reanalysis (JRA25; Onogi et al. 2007). Monthly mean river runoff is 
based on the Arctic Runoff Data Base (ARDB) as prepared by P. Winsor (2007, personal 
communication). 

The ECC02 Bering Strait solution was optimized based on observations (Nguyen 
et al. 2010). The model bathymetry in the vicinity of Bering Strait and the location of the 
Bering Strait cross-section are shown in Figure 29c. 

3. Naval Postgraduate School Arctic Modeling Effort (NAME) 

The NAME coupled sea-ice-ocean model was described earlier in Chapter II. 
There are 15 grid cells across Bering Strait in this model (Figure 29d). 

4. Nucleus for European Modelling of the Ocean (NEMO) With ORCA 
Configuration 

The ORCA025-N102 model configuration of the National Oceanography Centre 
Southampton is an “eddy-permitting” z-level global coupled sea ice-ocean model. 
ORCA025-N102 was developed within the Nucleus for European Modelling of the 
Ocean (NEMO) framework for ocean climate research and operational oceanography 
(http://www.nemo-ocean.eu/; Madec 2008) as part of the DRAKKAR configurations 
(DRAKKAR group 2007) and is largely based on the ORCA025-G70 configuration (e.g., 
Lique et al. 2009). ORCA025-N102 includes the ocean circulation model OPA9 (Madec 
et al. 1998) coupled to the Louvain-la-Neuve Ice Model sea ice model LIM2 (Fichefet 
and Morales Maqueda 1997). The ocean model is configured on a tri-polar Arakawa C- 
grid (Arakawa 1966) with the model poles at the geographical South Pole, in Siberia and 
in the Canadian Arctic Archipelago. The horizontal resolution is approximately 28 km at 
the equator, increasing to 6-12 km in zonal and ~3 km in meridional directions in the 
Arctic Ocean. The model resolves large eddies (-30-50 km), while “permitting” most of 
smaller eddies. ORCA025-N102 has a higher vertical resolution than the ORCA025-G70 
configuration, utilizing 64 vertical levels with thicknesses ranging from approximately 6 
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m near the surface to 204 m at 6000 m. The high vertical resolution in the upper ocean (8 
levels in the upper 50 m and 13 levels in the upper 100 m) greatly improves the model 
representation of the shallow Arctic continental shelves, Bering and Chukchi Seas. There 
are 8 model cells across Bering Strait (Figure 29e). The fine model resolution in the both, 
horizontal and vertical, together with high-resolution model bathymetry adapted from 
ETOP02 and partial steps in the model bottom topography accurately approximates the 
steep seabed relief near the Arctic shelves, resulting in the more realistic along-shelf flow 
(e.g., Bamier et al. 2006; Penduff et al. 2007). The LIM2 sea ice model uses the Viscous- 
Plastic (VP) ice rheology (Hibler 1979) and the 3-layer Semtner (1976) thennodynamics 
updated with sub-grid scale sea ice thickness distribution (Fichefet and Morales Maqueda 
1997) and sea ice thickness-dependent albedo (Payne, 1972). To obtain more distinct sea 
ice edges, the model employs the positive-definite, second moments conserving 
advection scheme by Prather (1986). The sea ice model is coupled to the ocean model 
every five oceanic time steps through a non-linear quadratic drag law (Timmermann et al. 
2005). 

For the 1958-2001 simulations used in the present study, the ORCA025 model 
was driven by the DRAKKAR Forcing Set (DFS 3.1) atmospheric reanalysis (Brodeau et 
al. 2010). The reanalysis combines monthly precipitation, daily downward shortwave and 
longwave radiation from the CORE forcing data set (Large and Yeager 2004) and 6- 
hourly 10 m wind, 2 m air humidity and 2 m air temperature from ERA40 reanalysis. The 
turbulent exchanges between atmosphere and ocean and atmosphere and sea ice are 
computed during model integration using the bulk formulae from Large and Yeager 
(2004). Climatological monthly continental runoff (Dai and Trenberth 2002) is included 
as an additional freshwater source, applied along the coastline. Initial conditions for 
temperature and salinity are derived from a monthly climatology that merges the Levitus 
(1998) World Ocean Atlas climatology with the PHC2.1 database (Steele et al. 2001) in 
high latitudes. To avoid salinity drift, the sea surface salinity is restored toward the 
monthly mean climatological values on the timescale of 180 days for the open ocean and 
12 days under sea ice. 
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5. Pan-Arctic Ice-Ocean Modeling and Assimilation System (PIOMAS) 

PIOMAS is a variant of BESTMAS (see description above) with a coarser 
horizontal resolution (~22 km) and smaller model domain (north of 49°N; Zhang et al. 
2008). However, it has 12 categories each for ice thickness, enthalpy, and snow depth 
(Zhang et al. 2008). The centers of the 12 ice thickness categories are 0, 0.26, 0.71, 1.46, 
2.61, 4.23, 6.39, 9.10, 12.39, 16.24, 20.62, and 25.49 m. The model bathymetry in the 
vicinity of Bering Strait and the location of the Bering Strait cross-section are shown in 
Figure 29f. 

C. BERING STRAIT OBSERVATIONAL MOORING DATA 

Year-round moorings have been deployed in the strait almost continuously since 
1990 (see Woodgate et al. 2006, 2010; and http://psc.apl.washington.edu/BeringStrait 
.html), at two locations, A2 and A3, as shown in Figure 29a. Site A2 is in the eastern 
portion of the strait and site A3 is just north of the strait (and in the US EEZ), and is 
hypothesized to provide a useful average of the flow through both of the channels 
(Woodgate et al. 2005a, b, 2006, 2007). For some years (92-93, 93-94, 94-95) the A3 
mooring was deployed ~ 120 nm further north, but these data are not considered here. 
Observations from A2 and A3 are available since autumn 1990, except for a few missing 
months, and for the deployment year autumn 1996-1997 when no moorings were 
deployed in the strait. In addition, measurements from site Al in the western channel of 
the strait, and thus in the Russian EEZ, are available for the early 1990s (from 90-91; 92- 
93, 93-94) and since 2004. A fourth mooring site, A4, was established near the U.S. 
coast in 2001 to measure the Alaska Coastal Current (Woodgate and Aagaard 2005). 
Measurements from these mooring sites were provided to the author for comparison with 
model results. A high-resolution array was deployed in the strait starting in 2007; for 
more details see http://psc.apl.washington.edu/BeringStrait.html. 

Since the region is ice-covered in winter, all mooring instrumentation has 
traditionally been kept near-bottom to avoid damage by ice keels. The moorings provide 
measurements of temperature, salinity and velocity at a single depth, approximately 10 m 
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above bottom. High correlation (0.95; Woodgate et al. 2005b) in velocity was found 
between all sites in the strait region (Woodgate et al. 2005b). Available Acoustic 
Doppler Current Profiler (ADCP) data (a few moorings and ship-based ADCP sections 
from the eastern channel), show strong coherence in the vertical (see e.g., Roach et al. 
1995), with some surface intensification of the flows, especially within the ACC. Based 
on spatial coherence, Woodgate et al. (2005a) used the near-bottom velocity 
measurements from A3 to estimate the total volume transport through Bering Strait. In 
tenns of water properties, clearly the near-bottom data do not capture the upper layer, 
which in the summer/autumn period of the year is likely 10-20 m thick, about 1-2 °C 
warmer and about 1 psu fresher than the lower layer (Woodgate and Aagaard 2005; 
Woodgate et al. 2010). 

The flow through the Bering Strait is generally believed to be driven by some far 
field forcing (often described as the pressure head forcing) modulated by local wind 
effects (see Woodgate et al. 2005b for discussion and historic references). Woodgate et 
al. (2005b) suggest this large-scale forcing likely explains the high correlation in 
velocities between sites. Near the Alaskan Coast on the edge of the eastern channel, 
there is a strong, surface-intensified seasonal current. This is the Alaskan Coastal 
Current, which is present from midsummer until about the end of the year (Paquette and 
Bourke 1974; Ahlnas and Garrison 1984; Woodgate and Aagaard 2005), and in summer 
CTD sections it is present as a ~10 km wide, 40 m deep wann, fresh current (Woodgate 
and Aagaard 2005). Much less is known about the Siberian Coastal Current, which is 
present sometimes on the Russian coast (Weingartner et al. 1999) and could influence the 
flow in the western part of the strait. 

D. RESULTS 

Model representations of the geographical width across Bering Strait range from 
90-160 km (Figures 30, 31). ORCA and PIOMAS have widths most similar to reality 
(—85 km), while BESTMAS, NAME, and ECC02 are wider than reality. In the model 
results (and somewhat in contradiction to observations), the long-term mean velocity 

structure across Bering Strait shows horizontal shear (Figure 31). The various horizontal 
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resolutions from the 5 models and the different bathymetry schemes make the results 
appear disparate upon first glance. In fact, the model cross-sectional area of the strait 
varies from 2.37-4.50 km , which is an ~8% smaller area in two models (NAME and 
PIOMAS) and 25-73% larger area in the other models compared to the observational 
area estimate of 2.6 km (Table 8). However, a closer look reveals agreement that 
horizontal shear is present and that the highest speeds tend to be in the eastern channel. 
Woodgate et al. 2005b showed no significant differences in the September 1990-October 
1991 observed mean near-bottom velocities between the two channels; however, at 
shorter time scales, the monthly mean velocities may differ by up to 35 cm/s. The multi¬ 
year mean flow from all models is northward, however the NAME model shows a weak 
(< 5 cm/s) reversal near the center of the strait. Note that Coachman et al. (1975), 
sometimes observe a minimum (and even reversal) in this region, although their sections 
could be subject to time aliasing as discussed below. Vertical shear is present in some 
model results, particularly the NAME, ECC02, and ORCA models. 
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Depth (m) Depth ( m ) 


BESTMAS (1979-2004) 



ORCA (1979-2001) PIOMAS (1979-2004) 



Figure 30. Vertical section of the long-term mean northward velocity (cm/s) across 
Bering Strait from all models. Positive velocity is northward. A black X 
marks the approximate location of the A2 mooring within each model 

domain. 
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Table 8. Cross-sectional area across Bering Strait for the models and observations 

and friction coefficients for the models. 




Bottom 

Lateral 

Lateral 

Surface Friction 
Coefficient 

Model/ 

Observed 

Area 

(km 2 ) 

Friction 

Coefficient 

Boundary 

Condition 

Friction 

Coefficient 

Ice- 

Ocean 

Air-Ocean 

BESTMAS 

3.24 

quadratic 
bottom drag: 
1.225 x 10' 3 

no-slip 

variable 

momentum 

harmonic 
horizontal 
mixing 
depending on 
variable grid 
size 

5.5 x 10' 3 

1.0 x 10' 3 

ECC02 

4.50 

quadratic 
bottom drag: 
2.1 x 10' 3 

free-slip 

modified Leith 
[Fox-Kemper 
and 

Menemenlis, 

20081 

5.4 x 10' 3 

Large and 
Pond 
(1981, 
1982) 

NAME 

2.37 

quadratic 
bottom drag: 
1.225 x 10' 3 

no-slip 

momentum 
biharmonic 
horizontal 
mixing: 
-1.25 x 10 18 
cm 4 /s 

5.5 x 10' 3 

0.6 x 10' 3 

Observed 

2.6 

N/A 

N/A 

N/A 

N/A 

N/A 

ORCA 

4.17 

quadratic 
bottom drag: 
1.0 x 10' 3 

free-slip 

bi-harmonic (- 
1.5e+ll m 4 /s) 

quadratic, 
5.0 x 10' 3 

CORE 

bulk 

formulae, 
Large and 
Yeager 
(2004) 

PIOMAS 

2.38 

quadratic 
bottom drag: 
1.225 x 10' 3 

no-slip 

variable 
momentum 
harmonic 
horizontal 
mixing 
depending on 
variable grid 
size 

5.5 x 10' 3 

1.0 x 10' 3 























































To compare with observations, monthly mean northward near-bottom velocity at 
sites A2 and A3 for models (color) and data (black) for 1979-2004 is shown in Figure 31. 
For A2 (eastern channel; Figure 29a), model velocities range from ~5 cm/s southward to 
over 80 cm/s northward. Predominantly, the flow is northward with the mean northward 
velocity ranging from 28.6 (+/- 1.0) to 40.1 (+/- 1.9) cm/s among models, over the time 
period when observations are available (Table 9). The range is 29.5 (+/- 0.49) to 43.2 
(+/- 0.88) cm/s over the larger 1979-2004 time period (Table 10) Two of the lower 
resolution models (ORCA and ECC02) have the highest velocities, while the higher 
resolution models (BESTMAS and NAME) have lower velocities. The observed mean 
northward velocity is 26.2 (+/- 2.8) cm/s, which is near the lower range of the modeled 
mean values. All of the models show a significant (at the 99% level) correlation with the 
observed velocities at this location. The correlation coefficients range from 0.67-0.78 for 
the monthly means at A2 (Table 11). 
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Velocity (cm/s) ry Velocity (cm/s) 


Velocity from ~1 Om above bottom at A2 location 


a) 



1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 


Velocity from ~10m above bottom at A3 location 



30.6 

29.9 

26.9 

20.9 

5.6 


BESTMAS - ECC02 - NAME - ORCA - PIOMAS-Data 

Figure 31. Monthly mean velocity at ~ 10 above the bottom from the A2 mooring 

location (upper) and A3 mooring location (lower). Model results are shown 
in color and the observations are shown in black. Mean values for the time 
period when data are available are shown on the far right. 
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Table 9. Mean velocity, volume transport, near-bottom temperature, and near-bottom 
salinity from the models and from observations, for the time period when 
observations are available, as shown in Figures 31-34. However, the values 
from EC002 are for 1992-2004 only and the values from ORCA are for 
1990-2001 only. Error estimates are shown in parenthesis. All model errors 
are calculated as the standard error of the mean (sample standard deviation 
divided by the square root of the sample size). 


Model/Data 

Mean 
velocity 
-10 m 
above 
bottom 
at the 
A2 

location 

(cm/s) 

Mean 
velocity 
-10 m 
above 
bottom 
at the 
A3 

location 

(cm/s) 

Mean 

volume 

transport 

(Sv) 

Mean 
temper¬ 
ature 
-10 m 
above 
bottom 
at the 
A2 

location 

(°C) 

Mean 
temper¬ 
ature 
-10 m 
above 
bottom 
at the 
A3 

location 

(°C) 

Mean 
salinity 
-10 m 
above 
bottom 
at the 
A2 

location 

(psu) 

Mean 
salinity 
-10 m 
above 
bottom 
at the 
A3 

location 

(psu) 

BESTMAS 

33.1 

(-1.3) 

■ 

m 

0.79 

(-0.25) 

0.65 

(-0.24) 


33.2 

(-0.06) 

Data 

26.2 

(-2.8) 

20.9 

(-2.3) 

0.8 

(0.20) + 

0.27 

(-0.30) 


32.26 

(-0.08) 


ECC02 

39.4 

(-1.2) 

26.9 

(-0.7) 

ESI 

1.1 

(-0.27) 

0.62 

(-0.23) 


32.29 

(-0.04) 

NAME 

33 

(-1.0) 

5.6 

(-0.5) 

m 

■ 

-1.26 

(-0.05) 

32.45 

(-0.04) 

32.61 

(-0.04) 

ORCA 




0.96 

(-0.25) 

0.79 

(-0.24) 

32.3 

(-0.05) 

32.47 

(-0.03) 

PIOMAS 

28.6 

(-1.0) 

29.9 

(-0.9) 

m 

0.34 

(-0.21) 

0.6 

(-0.23) 

M 

32.8 

(-0.05) 


+the uncertainty for the data estimate is -25% (Woodgate et al. 2005a, b) 
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Table 10. Long-term mean velocity, volume transport, near-bottom temperature, and 
near-bottom salinity from the models, for the time periods shown in Figures 
31-34. Error estimates are shown in parenthesis. 


Model/Data 

Mean 
velocity 
~10 m 
above 
bottom 
at the 
A2 

location 

(cm/s) 

Mean 
velocity 
~10 m 
above 
bottom 
at the 
A3 

location 

(cm/s) 

Mean 

volume 

transport 

(Sv) 

Mean 
temper¬ 
ature 
-10 m 
above 
bottom 
at the 
A2 

location 

(°C) 

Mean 
temper¬ 
ature 
-10 m 
above 
bottom 
at the 
A3 

location 

(°C) 

Mean 
salinity 
-10 m 
above 
bottom 
at the 
A2 

location 

(psu) 

Mean 
salinity 
-10 m 
above 
bottom 
at the 
A3 

location 

(psu) 

BESTMAS 

34 

22.8 



0.71 

33.18 


(-0.69) 

(-0.45) 


B!ftM 

(-0.16) 

(-0.05) 


ECC02 

39.9 

27.4 

mm 

1.08 

0.63 


32.25 

(-0.78) 

(-0.52) 


(-0.25) 

(-0.20) 


(-0.03) 

NAME 

34.1 

5.6 



-1.27 

32.46 

32.58 

(-0.52) 

(-0.24) 


(-0.06) 

(-0.03) 

(-0.02) 

(-0.02) 

ORCA 

43.2 

31.6 


0.89 

0.69 

32.33 

32.47 

(-0.88) 

(-0.65) 


(-0.17) 

(-0.15) 

(-0.02) 

(-0.02) 

PIOMAS 

29.5 

29.2 


0.26 

0.53 



(-0.49) 

(-0.49) 


(-0.14) 

(-0.15) 




Table 11. Correlation coefficients between models and the observations of northward 
velocity, temperature, and salinity at A2 and A3 locations. All correlations 
are significant at the 95% level. 


Model 

Velocity 

Temperature 

Salinity 

A2 

A3 

A2 

A3 

A2 

A3 

BESTMAS 

0.78 

0.71 

0.78 

0.7 

0.67 

0.53 

ECC02 

0.67 

0.72 

0.77 

0.76 

0.6 

0.48 

NAME 

0.69 

0.27 

0.73 

0.86 

0.7 

0.57 

ORCA 

0.68 

0.7 

0.79 

0.76 

0.6 

0.39 

PIOMAS 

0.70 

0.82 

0.88 

0.79 

0.66 

0.59 
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Figure 31b shows the near-bottom northward velocity at the A3 location. The 
model spread of velocities is larger for A3, with the NAME model having the lowest 
mean velocity (5.6 +/- 0.5 cm/s) and ORCA having the highest mean velocity (30.6 +/- 
1.3 cm/s) over the same time period as observations. The observed mean northward 
velocity is 20.9 +/- 2.3 cm/s. The NAME model has a problem representing the correct 
magnitude of the velocity at this location, which appears to be related to the upstream 
bathymetry representation. The correlation between the models and the data is more 
variable at A3, with correlation coefficients ranging between 0.27 and 0.82 (Table 11). 

Near-bottom monthly mean temperatures at the A2 and A3 mooring locations are 
shown in Figure 32. (Temperature at A1 is not shown because there are too few data 
available at this time.) Temperatures tend to be warmer at the southern A2 location, with 
model means ranging between -0.96 (+/- 0.09) to 1.1 (+/- 0.27) °C. The mean observed 
near-bottom temperature for the same location is 0.27 (+/- 0.3) °C. ORCA, ECC02, and 
BESTMAS models tend to overestimate the temperature by 0.5-0.8 °C in the mean, while 
NAME underestimates the temperature by 1.2 °C in the mean. It is possible that the 
colder temperatures for the NAME model may be related to excessive ice production, 
which causes the sinking of very cold surface water, especially in polynya regions of the 
northern Bering Sea. Surprisingly, the PIOMAS temperatures are closest to the observed, 
despite the fact that it is the lowest resolution model in this study and only has 3 grid 
points across the strait (Figure 30). Temperatures at the A3 location are, again, 
underestimated in the NAME model and overestimated in ORCA, ECC02, BESTMAS 
and also in PIOMAS. While the magnitude of the model-data differences may be up to 
~1 °C in the mean, the models’ results are significantly correlated (at the 99% confidence 
level) with the observations. The correlation coefficients range between 0.73-0.88 at A2 
and between 0.70-0.86 at A3 (Table 11). There is no trend, either observed or modeled, 
in the time series shown here. There is, however, a strong annual cycle present (as 
discussed by Woodgate et al. 2005a), which dominates the correlations. This annual 
cycle will be discussed below. 
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Mean Temperature at ~10m above bottom at A2 location 




Mean Temperature at ~10m above bottom at A3 location 



BESTMAS - ECC02 - NAME - ORCA - PIOMAS-Data 


Figure 32. Monthly mean near-bottom temperature (°C) at the (a) A2 and (b) A3 
mooring locations. Model results are shown in various colors and 
observations are shown in black. Mean values for the time period when data 
are available are shown on the far right. 


A similar analysis was perfonned for salinity at the A2 and A3 mooring locations 
(Figure 33). The mean modeled salinity ranges between 31.7 (+/- 0.06) and 33.2 (+/- 
0.06) psu at A2 and between 32.2 (+/- 0.04) and 33.2 (+/- 0.06) psu at A3. The mean 
observed salinities are 32.3 (+/- 0.08) at A2 and 32.5 (+/- 0.06) at A3. The BESTMAS 
and PIOMAS models tend to overestimate the salinity, by up to 0.9 psu above the 
observed mean value, whereas the NAME, ECC02, and ORCA models have values close 
to the observed. All of the models’ results are significantly correlated (at the 99% 
confidence level) with the observations of salinity at A2 and A3. The correlation 
coefficients range between 0.60 and 0.70 at A2 and between 0.39 and 0.59 at A3 (Table 
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11). The correlations are not as high for salinity as they are for temperature, especially at 


the A3 location. Again, an annual cycle of salinity is apparent in the time series (also see 
Woodgate et al. 2005a); however, it is not as strong as the annual cycle of temperature. 
The actual annual cycles will be discussed below. 


a) 


Mean Salinity at ~10m above bottom at A2 location 



33.09 

32.45 

32.30 

32.26 

31.72 


Mean Salinity at ~10m above bottom at A3 location 



BESTMAS - ECC02 - NAME - ORCA - PIOMAS-Data 


Figure 33. Monthly mean near-bottom salinity at the (a) A2 and (b) A3 mooring 
locations. Model results are shown in various colors and observations are 
shown in black. Mean values for the time period when data are available are 

shown on the far right. 


It is important to recognize that a comparison between point measurements and 
model results is difficult. In the data, velocity, as well as temperature and salinity, is 
measured at a single point, while in models it is a grid-cell mean, which may range from 
a few to tens of kilometers in the horizontal and several meters in the vertical. In addition, 
the discrepancy between the real and model bathymetry introduces a difference in 
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bathymetric gradients, displacing model currents from their “true” geographical 
positions. Table 12 and Figure 30 illustrate these points. Table 12 shows the depth at the 
moorings A2 and A3 and model depth at the co-located virtual moorings; the difference 
between real and model bathymetry is clear. Moreover, in the models (although 
seemingly less so in the observations) velocity can vary significantly between the 
adjacent model grid cells (Figure 30), thus the results of model-data comparisons would 
strongly depend upon the exact geographical position of model virtual moorings. Finally, 
the stochastic nature of the oceanic turbulence is not simulated by the models used in this 
study. Therefore, it is likely more informative to evaluate model results using integrated 
fluxes, as discussed below. 

Table 12. Depth information (m) for the models and the observations at the A2 and 

A3 mooring locations. 


Location 

Model/data 

Water 

column 

depth 

(m) 

Mid-depth of 
model grid cell 
or depth of 
observation ~10 
m above 
bottom (m) 

A2 

Data 

53.0 

44.0 

BESTMAS 

51.0 

39.5 

ECC02 

50.0 

35.0 

NAME 

53.0 

37.7 

ORCA 

57.9 

35.5 

PIOMAS 

43.0 

33.0 

A3 

Data 

56.0 

47.0 

BESTMAS 

51.0 

39.5 

ECC02 

50.0 

35.0 

NAME 

53.0 

37.7 

ORCA 

57.9 

35.5 

PIOMAS 

43.0 

33.0 
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Monthly mean Bering Strait volume transport from the models and observations 
is shown in Figure 32. Model means range from 0.67 (+/- 0.03) to 1.29 (+/- 0.06) Sv 
(Table 9) over the time period when observations are available. The volume transport is 
highest for the ORCA and ECC02 models and is lowest for the PIOMAS, BESTMAS 
and NAME models. The observed estimate of the long-term mean (1991-2004) volume 
transport through Bering Strait is 0.8 +/- 0.2 Sv (Woodgate et al. 2005a). This estimate is 
based on single-depth velocities at the A3 mooring location multiplied by the total cross- 
sectional area of 4.25 km", although numbers do not differ significantly if using 
observations from the A2 site. 

Bering Strait Volume Transport 


t— i—i—i—i—i—i—i—i—n—i—i—i—i—i—i—i—i—i—i—r 



1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 


BESTMAS - ECC02 - NAME - ORCA - PIOMAS-Data 

Figure 34. Monthly mean volume transport from the models and observations. The 
observations are based on the near-bottom velocity at the A3 mooring 
location multiplied by a cross-sectional area of 4.25 km", as per Woodgate 
et al. (2010). Mean values for the time period when data are available are 

shown on the far right. 


Annual mean volume transport from models and observations is shown in Figure 
35a. Observed volume transport ranges from 0.6-1 Sv (+/- 0.2 Sv; Woodgate et al. 
2006), which is most similar to the estimates from the BESTMAS, NAME, and PIOMAS 
models. The ACC, which is not observed by the A2 or A3 near-bottom measurements, 
may add around 0.1 Sv to the summer estimates (Woodgate and Aagaard 2005); 
however, the main uncertainty with these estimates is the assumption of spatially uniform 
flow across the strait. 


98 












































Bering Strait Volume Transport Annual Means 




Bering Strait Freshwater Flux Annual Means 



BESTMAS ECC02 — *— NAME ORCA — * — PIOMAS 

-»— Data (without ACC estimate) .«•.Data (with ACC estimate) 


Figure 35. Annual mean (a) volume transport, (b) heat, and (c) freshwater fluxes. Heat 
is referenced to -1.9 °C for NAME, ECC02, and ORCA, in order to 
compare with cited observations in the text. Heat is referenced to -1.8 °C 
for PIOMAS and BESTMAS, due to model limitations. Freshwater is 
referenced to 34.8 psu. Observed volume transport and freshwater flux 
values do not include the ACC and stratification, which likely add ~ 0.1 Sv 
and 800-1,000 km 3 /yr (see Woodgate et al. 2006). The dashed black line 
(c) represents the freshwater flux with the estimated ACC and stratification 
correction of an additional 900 km 3 /yr. The observed heat flux values 
include an estimate for the ACC using SST for a 10-m surface layer (lower 
bound; lower black line) and a 20-m surface layer (upper bound; upper 
black line). Observed heat flux values are described further in Woodgate et 

al. (2010). 
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Freshwater flux from the Bering Sea into the Chukchi Sea is an important factor 
affecting stratification and the maintenance of the Arctic Ocean halocline (e.g., Aagaard 
et al. 1985a). Heat fluxes through Bering Strait and through the Chukchi shelf appear to 
influence the distribution and thickness of sea ice (Shimada et al. 2006; Woodgate et al. 
2010). Previously published observations of heat flux (e.g., Woodgate et al. 2010) use a 
reference temperature of -1.9°C, which will overestimate heat fluxes when freezing 
temperature is above that value. However, for comparison, the same value was used for a 
reference temperature in the model calculations. The PIOMAS and BESTMAS models 
use constant reference temperature of -1.8°C to calculate heat flux because these two 
models use this temperature as the freezing temperature for an ease in conserving heat in 
the model. It is estimated that this difference in reference temperature between -1.9°C 
and -1.8°C produces errors that are < 5% of the mean. Oceanic heat flux through Bering 
Strait in the models was calculated as the vertical and horizontal integral of: the heat 
(heat capacity multiplied by the difference between the temperature and the reference 
temperature) multiplied by grid cell velocity normal to the cross-section on a monthly 
mean time scale. 

The annual mean oceanic heat flux time series for the models and observations (as 
per Woodgate et al. 2010) are shown in Figure 37b. In the models, peaks in the heat flux 
occurred during several years (e.g., 1979, 1986, 1993, and 1997) and consistently showed 
up in results from all five models. However, data coverage is not sufficient to confirm 
these peaks. A peak in 2004 is noted in observations (see Woodgate et al. 2010) and is 
apparent in all of the models, except ORCA, which does not have results for that time 
period. ECC02 is also able to simulate a recent increase in heat flux in 2007 (not shown), 
similar to the observations (Woodgate et al. 2010). 

20 

The long-tenn model mean heat flux ranged between 1.5-5.1 x 10 J/yr. ORCA 
and ECC02 have much higher values than BESTMAS and NAME. Observations of the 
annual heat flux based on near-bottom measurements, a correction for the ACC, and SST 
from satellite data were published in Woodgate et al. (2010). The observationally 
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estimated range of heat flux is ~2.8-4.5 x 10 20 J/yr with estimated uncertainty of 0.8 x 
10 20 J/yr, based on years 1991, 1998, 2000-2006. However, the 2007 heat flux was 
estimated at over 5x10“ J/yr. 

For the calculation of freshwater fluxes, a reference salinity of 34.8 psu was used 
because this value is considered to be the mean salinity of the Arctic Ocean and has been 
used in most other Arctic studies (based on original work by Aagaard and Carmack 
1989). Integrated annual mean oceanic freshwater fluxes were calculated on a monthly 
mean timescale (see Equation 1 in Melling 2000) from each of the models and are shown 
in Figure 37c. An observationally-based lower bound of annual mean freshwater fluxes 
is also shown, however these values do not include the ACC or stratification and 
according to Woodgate et al. (2005c, 2006) they underestimate the freshwater flux by 
about 800-1,000 km /yr. With this correction, the observed freshwater annual means are 
similar to results from the ECC02 and ORCA models, with the other models appearing 
to underestimate the total freshwater flux. No long-term trend is apparent in either the 
heat or freshwater flux for this time period, however a gradual increase in freshwater 
during the early 2000s has occurred in the model results, ending with a peak in freshwater 
flux in 2004, similar to observations (also see Woodgate et al. 2006, 2010). 

It is important to note that both the models and the data have limitations with 
respect to calculations of heat and freshwater fluxes. The models used here are too 
coarse to represent the narrow Alaska Coastal Current (ACC), which is estimated to carry 
10-20% of the freshwater flux (Woodgate and Aagaard 2005) and 20-30% of the heat 
flux (Woodgate et al. 2010) through the strait. The historic near-bottom data used here 
does not measure the upper ocean seasonally varying stratification or ACC, which is a 
surface/coastal feature. Thus, on-going observational research is using extra moorings, 
hydrographic data and upper water column sensors to estimate stratification (see e.g., 
http://psc.apl.washington.edu/BeringStrait.html). Also, new models are configured at 
increasingly high spatial resolution, an example of which is presented at the end of this 
chapter. 


101 



Arctic shelf seas have a strong annual cycle of temperature and salinity, especially 
in the upper ocean; some areas may also exhibit strong annual changes in the oceanic 
circulation. The Bering Strait region is no different in this respect. Observations have 
shown stronger northward flows in summer (e.g., Woodgate et al. 2005a, b and 
references therein). According to the model results, volume transport peaks in summer 
(May-July) and is lowest in winter (December-March; Figure 36) and this agrees 
reasonably with observational results (peaking in May/Jun, minimum in December- 
February; Woodgate et al. 2005a). However, error bars make these maxima less certain in 
both the data and model results. The point data measuring near-bottom flow have a 
larger annual cycle, with a range of 0.4 to 1.3 Sv (errors order 25%; Woodgate et al. 
2005a, b). PIOMAS, BESTMAS, ECC02, and NAME models combining all points 
across the strait have similar annual cycles to the data; however, they are not as strong. 
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a) 


Bering Strait Volume Transport Annual Cycle (1991 -2004) 





BESTMAS —*— ECC02 —•— NAME —•— ORCA —•— PIOMAS —Data 


Figure 36. Annual cycles of (a) volume transport, (b) heat flux, and (c) freshwater 
transport. The annual cycles are averaged over 1991-2004, except for 
ORCA (1991-2001) and ECC02 (1992-2004). The heat flux is referenced 
to -1.9 °C for NAME, ECC02, and ORCA. The reference temperature is - 
1.8°C for PIOMAS and BESTMAS due to model limitations. The 
freshwater transport is referenced to 34.8 psu. 


103 






























































As shown by Woodgate et al. (2010), the heat flux annual cycle is also very 
strong. Observational results from near-bottom moorings combined with SST from 
satellite data (Woodgate et al. 2010, Figure 3) suggest strong interannual variability in the 
timing of the summer peak, although the computation presented there does not include 
the seasonality of the ACC. In the models (Figure 36b), heat flux peaks in summer and is 
near-zero in winter. However, the models do not agree on the magnitude of the 
summertime peak, which ranges between 15 (+/- 6.4) to over 40 (+/- 14) TW. The heat 
flux is near zero for December-April (when water temperatures are around freezing). 
The models with the highest resolutions (BESTMAS and NAME) show lower peaks in 
the summertime heat flux [15 (+/- 6.4) and 22.5 (+/- 7.9) TW], while the lower resolution 
models have higher heat fluxes, which is in part related to their simulation of summer 
volume flux. 

Freshwater flux is important to the Chukchi Sea and the Arctic Ocean due to the 
importance of maintaining the halocline there. This halocline separates the upper mixed 
layer from the deeper, warmer Atlantic layer. Annual cycles of freshwater flux through 
Bering Strait are similar for PIOMAS, BESTMAS, and NAME, with peaks in the 
summer (June-August) and lowest in winter (December-April; Figure 36c). Again, error 
bars make these peaks in the model results less certain. The freshwater flux maxima for 
these models is between 65 (+/- 14.3) - 80 (+/- 13.4) mSv in July. Annual cycles for 
ECC02 and ORCA have somewhat similar shapes, however they transport more 
freshwater (up to 115 (+/- 11.7) mSv in summer and more than 60 (+/- 40) mSv in winter 
for ORCA) to the north. 

Using the monthly mean volume transport values shown in Figure 36 and the 
annual cycles of volume transport shown in Figure 36a, I calculated the volume transport 
anomalies, which are shown in Figure 37. The anomalies confirm that there is no trend in 
the transport time series from 1979-2004. The ORCA model tends to produce larger 
anomalies than any of the other models or data. However, all of the models are 
significantly correlated (at the 95% level) with the data. The correlation coefficients 
range from 0.59 to 0.73. 
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Bering Strait Volume Transport Anomaly 
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Figure 37. Monthly mean volume transport anomalies from the annual cycle (Sv). The 
annual cycle for each model and the data are shown in Figure 36a. 

E. SUMMARY AND CONCLUSIONS 

Model volume transports ranged from 0.67 (+/- 0.03) to 1.29 (+/- 0.06) Sv in the 
mean, compared to observational estimates of 0.8 +/-0.2 Sv. Thus, most of the models 
are in agreement with the observational estimate to within errors. ORCA and ECC02 
showed the highest volume transports, while NAME and BESTMAS showed the lowest 
transports. Note that of the models, ORCA and ECC02 also have the largest cross 
sectional area of the strait. Higher resolution models seem to give lower transport 
estimates; the cause of this is not fully understood. The models are using both lateral and 
vertical friction parameterization to represent the flow next to the boundary 
(bottom/surface or lateral; see Table 8). Some uncertainty of model estimates of the 
volume transport throughout the strait might be related to the estimation of the frictional 
layers, subject to the parameterization used. Using a 1/4° global model, Penduff et al. 
(2007) have demonstrated that enstrophy-conserving momentum advection schemes 
produce a spurious numerical sidewall friction, leading to a weaker topographic 
alignment of the mean flow and weaker barotropic transports. Both spurious sidewall 
and explicit no-slip lateral friction could contribute to lower transports in BESTMAS and 
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NAME models compared to ECC02 and ORCA, as the last two models feature free-slip 
lateral boundary conditions. Besides, ORCA utilizes energy-enstrophy conserving 
advection, resulting in higher transport than in ECC02 (Tables 9 and 10). However, this 
cannot explain a higher transport in the PIOMAS model compared to BESTMAS, since 
these two models share the same configuration, except for the resolution and different 
number of sea ice categories (12 and 8 respectively). 

Panteleev et al. (2010) applied an inverse model (with 10 grid points across the 
strait) to reconstruct the flow using available data for 1990-1991 and recently calculated 
the transport through Bering Strait as 0.57 Sv (no stated uncertainty). The data used to 
reconstruct the circulation were from 12 moorings that were deployed in the Bering Strait 
and Chukchi Sea from September 1990 to October 1991 (Woodgate et al. 2005b). This 
estimate from Panteleev et al. (2010) tends to agree with the estimates from the 
BESTMAS and NAME models. In fact, the mean volume transports during the same 
time period (September 1990-October 1991) were 0.62 (+/- 0.03) and 0.59 (+/- 0.03) Sv 
for the BESTMAS and NAME models, respectively. 

The model sections presented here show significant vertical and horizontal 
velocity shear across the strait. This appears to be somewhat in contrast to the limited 
observations, which show strong coherence of flow and agreement of near-bottom speeds 
in the 2 channels of the strait. However, a strong coherence between 2 point 
measurements in the strait does not indicate that the flow is homogenous throughout the 
strait. Observations indicate a stronger flow in the ACC (not shown), which would 
suggest some horizontal velocity shear across the strait. The only currently published 
sections of observed velocity in the strait are those of Coachman et al. (1975), which 
indicate both horizontal and vertical shear in the speeds across the strait. However, these 
sections may be subject to time aliasing being taken over a period of days. 

It seems possible that some of the horizontal and vertical shear found in the 
models is due to edge effects and/or the poor resolution of the real world bathymetry. It 
must also be remembered however, that the observational transports presented here are 
based on near-bottom measurements of velocity at A2 or A3 and an assumption of 
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homogeneity of flow at all locations in the strait, which may not be an accurate 
assumption. This is currently being tested by an increased mooring effort in the strait 
region. Future observational results from this expanded mooring array should help to 
narrow the uncertainty in the measurements and produce more reliable integrated fluxes. 

Increases in model resolution and improvements in model bathymetry are 
expected to improve the representation of the flow structure in Bering Strait. Initial spin- 
up results from a higher resolution version (-2.3 km in the horizontal) of the NAME 
model are compared with the current 9-km results in Figure 37. The improvements in the 
strait bathymetry are obvious, with actual Diomede Islands being present in the strait and 
a more realistic depth of -50 m is reached in the two channels of the strait. However, 
representation of the northward velocity and its horizontal shear is not qualitatively 
different between the two models. The new 2.3-km model has 4 times as many (i.e., 60) 
grid cells across the strait, which allows for the representation of the narrow Alaska 
Coastal Current and its velocity core in the eastern part of the strait appears to be 
separated from the coast by a narrow boundary layer, which may be more similar to 
recent observations (Rebecca Woodgate, 2010, personal communication). When new 
model output becomes available, a future manuscript is planned to analyze those results 
in comparison with data from the extended array of moorings across the Bering Strait. 
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Figure 38. Bering Strait velocity (cm/s) from the (a) 9-km model during the 26-year 
mean March and (b) 2-km model during spin-up mean March. 


108 


















VIII. LINKAGES TO WEATHER INDICES, BIOLOGICAL 
PRODUCTIVITY, AND THE WESTERN ARCTIC 


A. LINKAGES TO LARGE-SCALE WEATHER INDICES 

Three large-scale weather indices of the Arctic and/or North Pacific have been 
li nk ed to several low-frequency oscillations within the northern North Pacific and the 
Bering Sea (e.g., Niebauer 1988; Niebauer and Day 1989; Niebauer et al. 1999). These 
indices are the Arctic Oscillation (AO), the Pacific Decadal Oscillation (PDO) and the 
Pacific/North American Pattern (PNA). 

The AO is an atmospheric circulation pattern in which the atmospheric pressure 
over the polar regions varies in opposition with that over middle latitudes (about 45° 
North) on time scales ranging from weeks to decades (see e.g., 
http://nsidc.org/arcticmet/pattems/arctic oscillation.html) . High pressure over the polar 
region and low pressure at mid-latitudes characterizes the “negative phase”, while the 
“positive phase” exhibits the opposite pattern. A positive AO tends to drive ocean storms 
farther north due to the higher pressure at mid-latitudes. 

The PDO Index is defined as the first principal component of extratropical North 
Pacific Ocean SST anomalies (Mantua et al. 1997; Zhang et al. 1997; Hare and Mantua 
2000). The PDO index is available at http://iisao.washington.edu/pdo/ . The PDO term 
was first introduced by Steven Hare in the late 1990s, who used it to help explain 
connections between Alaska salmon production cycles and Pacific climate (Hare et al. 
1999). Major changes in northeast Pacific marine ecosystems have been correlated with 
phase changes in the PDO (e.g., Beamish 1993; Hare and Mantua 2000). Warm 
(positive) time periods have seen enhanced coastal ocean biological productivity in 
Alaska and inhibited productivity off the west coast of the contiguous United States, 
while cold (negative) PDO time periods have been associated with an opposite north- 
south pattern of marine ecosystem productivity. 

The Pacific/North American teleconnection pattern (PNA) is an important mode 
of low-frequency variability in the Northern Hemisphere extratropics 
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( http://www.cpc.noaa.gov/data/teledoc/pna.shtnil ; Barnston and Livezey 1987). The 
positive phase of the PNA pattern exhibits above-average atmospheric geopotential 
height fields in the vicinity of Hawaii and over the Rocky Mountain region of North 
America, and below-average heights located south of the Aleutian Islands and over the 
southeastern United States (Wallace and Gutzler 1981). The negative phase is associated 
with a westward retraction of the jet stream toward eastern Asia, blocking activity over 
the high latitudes of the North pacific, and a strong split-flow configuration over the 
central North Pacific. Anomalously warm temperatures over western Canada and the 
western United States are associated with a positive phase of the PNA, as well as high 
levels of precipitation from the Gulf of Alaska into the Pacific Northwestern United 
States. Yin (1994) demonstrated that the PNA is also strongly influenced by the El Nino/ 
Southern Oscillation (ENSO) phenomenon. The positive phase of the PNA pattern tends 
to be associated with Pacific warm episodes (El Nino), and the negative phase tends to be 
associated with Pacific cold episodes (La Nina). 

The monthly mean and 13-month running mean values of the AO, PDO, and PNA 
are shown in Figure 39 for 1979-2004. A 13-month running mean was chosen in order 
to smooth out the annual cycle. 
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Climate Indices 



1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 2004 


Climate Indices (13-mo. running mean) 



AO - PDO - PNA 


Figure 39. Monthly mean (upper) and (lower) 13-month running mean values of the 

AO, PDO, and PNA. 

The low-frequency temperature changes in the Bering Slope Current (discussed in 
Chapter V) may be indirectly related to large-scale weather indices of the Arctic and/or 
North Pacific. Correlations between the monthly mean temperatures at sections BSC1, 
BSC2, BSC3, and BSC4 and the AO, PDO, and PNA were calculated for a lag range of 
+/- 50 months. A lag range of +/- 50 months was used because the low-frequency 
oscillations in Figure 39 are of a similar time period. Although the monthly mean 
temperatures are not significantly correlated at any of these BSC sections, the 13-month 
running means showed some significant correlations (Table 13). 
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Table 13. Correlation coefficients between the mean temperature across BSC sections 
and large-scale weather indices. Correlations are for the 13-month running 
mean values of both the temperature time series and the weather indices for 
1979-2004. Correlations with no lag are shown (Corr. no lag) with the 
associated p-value. The maximum correlation (Max. Corr.) is also shown 
for a lag range of +/- 50 months, along with the lag time and p-value. A 
negative lag means that the mean temperature across the BSC section is 
leading the weather index. The significance level used is 95%. 


Section 

Index 

Corr. no lag 

P-value 

Max. Corr. 

Lag 

P-value 

BSC1 

AO 

-0.33 

0.0000 

-0.57 

-8 

0.0000 

BSC1 

PDO 

0.02 

0.6764 

0.33 

50 

0.0000 

BSC1 

PNA 

0.15 

0.0080 

0.32 

-9 

0.0000 

BSC2 

AO 

-0.27 

0.0000 

-0.59 

-36 

0.0000 

BSC2 

PDO 

0.09 

0.1348 

0.38 

50 

0.0000 

BSC2 

PNA 

0.10 

0.0936 

0.42 

-14 

0.0000 

BSC3 

AO 

-0.29 

0.0000 

-0.52 

-17 

0.0000 

BSC3 

PDO 

0.09 

0.1178 

0.40 

-44 

0.0000 

BSC3 

PNA 

0.28 

0.0000 

0.30 

3 

0.0000 

BSC4 

AO 

-0.27 

0.0000 

-0.51 

-34 

0.0000 

BSC4 

PDO 

0.13 

0.0261 

0.32 

50 

0.0000 

BSC4 

PNA 

0.12 

0.0310 

0.30 

-50 

0.0000 


The AO is negatively correlated with the mean temperature across all of the BSC 
sections. Although the correlation coefficients are somewhat low (-0.27 to -0.33), they 
are statistically significant at the 95% level with a zero lag time. This suggests that a 
positive AO, which is associated with increased storminess or weakened high pressure in 
the central Arctic leads to lower temperatures along the slope in the Bering Sea; however 
such a connection, if it exists at all, would require further verification. Using a range of 
+/- 50 months lag time, some correlations are higher (up to -0.59); however, the 
interpretation of the lagged correlations is even less conclusive. Correlations of 
temperature across the BSC with the PDO/PNA do not appear to be consistent or 
sensible. 
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Correlations between the 13-month running mean net volume transport through 
Bering Sea straits/passes and large-scale weather indices (AO, PDO, and PNA) are 
shown in Table 14. Although there are some significant correlations between Aleutian 
Island passes and straits, a clear pattern does not emerge with any of the indices. The 
correlations are not consistent for a given weather index, instead they may be negative or 
positive along the Aleutian Island sections. This suggests that the inflow into the Bering 
Sea is not strongly related to large-scale weather patterns. As shown in Chapter IV, this 
inflow appears to be more closely related to the strength and position of the Alaskan 
Stream and mesocale eddies within the stream and in the southern Bering Sea. 
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Table 14. Correlation coefficients between the net volume transport through Bering 
Sea cross-sections and large-scale weather indices. Correlations are for the 
13-month running mean values of both the volume transport time series and 
the weather indices for 1979-2004. Correlations with no lag are shown 
(Corr. no lag) with the associated p-value. The maximum correlation (Max. 

Corr.) is also shown for a lag range of +/- 50 months, along with the lag 
time and p-value. A negative lag means that the volume transport through 
the cross-section is leading the weather index. 


Section 

Index 

Corr. no 
lag 

P-value 

Max. 

Corr. 

Lag 

P-value 

Unimak 

AO 

-0.13 

0.0201 

-0.35 

11 

0.0000 

Unimak 

PDO 

0.15 

0.0073 

-0.45 

-19 

0.0000 

Unimak 

PNA 

0.19 

0.0010 

0.41 

-44 

0.0000 

Akutan 

AO 

-0.08 

0.1765 

0.40 

-29 

0.0000 

Akutan 

PDO 

0.29 

0.0000 

0.56 

9 

0.0000 

Akutan 

PNA 

0.09 

0.1107 

0.42 

-46 

0.0000 

Urnnak 

AO 

-0.05 

0.3855 

0.39 

-34 

0.0000 

Urnnak 

PDO 

0.27 

0.0000 

0.50 

8 

0.0000 

Urnnak 

PNA 

0.07 

0.2027 

0.35 

-44 

0.0000 

StraitH 

AO 

-0.05 

0.3852 

0.18 

-43 

0.0046 

StraitH 

PDO 

0.05 

0.3836 

-0.38 

-19 

0.0000 

StraitH 

PNA 

0.07 

0.2227 

-0.20 

-23 

0.0006 

Samalga 

AO 

-0.03 

0.5671 

0.46 

-35 

0.0000 

Samalga 

PDO 

0.29 


0.45 

7 

0.0000 

Samalga 

PNA 

0.05 


-0.34 

-27 

0.0000 

StraitG 

AO 

-0.16 

0.0048 

-0.42 

44 

0.0000 

StraitG 

PDO 

0.23 

0.0001 

-0.37 

-21 

0.0000 

StraitG 

PNA 

0.14 

0.0168 

-0.34 

-29 

0.0000 

Amukta 

AO 

-0.27 


-0.45 

-33 

0.0000 

Amukta 

PDO 

0.17 


0.41 

-7 

0.0000 

Amukta 

PNA 

0.36 

0.0000 

0.43 

-7 

0.0000 

Seguam 

AO 

0.02 

0.7751 

0.46 

-34 

0.0000 

Seguam 

PDO 

0.24 

0.0000 

0.52 

8 

0.0000 

Seguam 

PNA 

0.06 

0.2690 

0.48 

-48 

0.0000 

StraitF 

AO 

-0.09 

0.1300 

0.36 

-35 

0.0000 

StraitF 

PDO 

0.27 

0.0000 

0.52 

8 

0.0000 

StraitF 

PNA 

-0.01 

0.8523 

0.37 

26 

0.0000 

StraitE 

AO 

-0.07 

0.2518 

0.43 

-35 

0.0000 

StraitE 

PDO 

0.39 

0.0000 

0.50 

5 

0.0000 
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Section 

Index 

Corr. no 
lag 

P-value 

Max. 

Corr. 

Lag 

P-value 

StraitE 

PNA 

0.01 

0.8479 

0.40 

26 

0.0000 

StraitD 

AO 

-0.05 

0.3979 

0.36 

-35 

0.0000 

StraitD 

PDO 

0.40 

0.0000 

0.48 

4 

0.0000 

StraitD 

PNA 

-0.01 

0.8926 

0.38 

26 

0.0000 

Amchitka 

AO 

0.02 

0.6677 

0.34 

49 

0.0000 

Amchitka 

PDO 

-0.35 

0.0000 

0.54 

-22 

0.0000 

Amchitka 


-0.09 

0.1097 

0.32 

38 

0.0000 

StraitC 

AO 

-0.13 

0.0287 

0.23 

-50 

0.0003 



0.18 

0.0017 

-0.52 

45 

0.0000 

StraitC 


-0.19 

0.0007 

0.35 

16 

0.0000 

Buldir 

AO 

-0.21 

0.0002 

-0.61 

11 

0.0000 

Buldir 

PDO 

0.22 

0.0002 

0.59 

10 

0.0000 

Buldir 


0.08 

0.1442 

-0.51 

-13 

0.0000 

StraitB 

AO 

-0.05 

0.4003 

-0.34 

22 

0.0000 

StraitB 

PDO 

-0.28 

0.0000 

0.38 

18 

0.0000 

StraitB 


-0.48 

0.0000 

-0.48 

0 

0.0000 

Near 1 

AO 

0.00 

0.9351 

0.39 

-27 

0.0000 

Near 1 

PDO 

0.08 

0.1718 

0.42 

29 

0.0000 

Near 1 

PNA 

-0.20 

0.0006 

-0.22 

3 

0.0001 

Near 2 

AO 

-0.16 

0.0054 

-0.30 

-35 

0.0000 

Near 2 

PDO 

0.13 

0.0212 

0.46 

-43 

0.0000 

Near 2 

PNA 

0.20 

0.0005 

0.32 

-46 

0.0000 

Near 

AO 

-0.15 

0.0112 

0.32 

19 

0.0000 

Near 

PDO 

0.16 

0.0050 

0.56 

-44 

0.0000 

Near 

PNA 

0.07 

0.1983 

0.31 

-48 

0.0000 

Kamchatka 

AO 

0.27 

0.0000 

0.27 

1 

0.0000 

Kamchatka 

PDO 

-0.11 

0.0485 

-0.47 

-43 

0.0000 

Kamchatka 

PNA 

-0.09 

0.1343 

-0.38 

-48 

0.0000 


AO 

-0.40 

0.0000 

-0.44 

-11 

0.0000 



0.02 

0.7045 

0.23 

12 

0.0001 



0.14 

0.0128 

0.40 

-43 

0.0000 


AO 

-0.17 

0.0028 

0.28 

-11 

0.0000 



0.02 

0.7862 

-0.35 

14 

0.0000 



0.22 

0.0001 

-0.31 

-50 

0.0000 


AO 

-0.49 

0.0000 

-0.51 

2 

0.0000 

Bering 

PDO 

0.03 

0.5978 

0.29 

-38 

0.0000 

Bering 

PNA 

0.28 

0.0000 

0.34 

-39 

0.0000 
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In the northern Bering Sea, there appears to be a consistent pattern of correlation 
with the AO. The net volume transports across Anadyr, Shpanberg, and Bering straits all 
show significant (at the 95% level) negative correlations with the AO at a zero time lag. 
The correlation is highest for Bering Strait (-0.49), with lower values at Anadyr (-0.40) 
and Shpanberg (-0.17). The correlations can be slightly higher when lag times are used, 
however they are not consistent among the sections, which minimizes their physical 
meaning. There is also an indication of a weak, but significant (at the 95% level), 
positive correlation between the northern straits and the PNA; correlation coefficients 
range between 0.14 and 0.28, with Bering Strait showing the highest correlation. 

B. PHYSICAL IMPACTS ON BIOLOGICAL PROCESSES 

The shelf waters of the northern Bering Sea are supplied with nutrients via 
upwelling of water from the deep Bering Sea basin (Tsunogai et al. 1979). The model 
results in Chapter IV provide a mechanism for how and where this upwelled water 
reaches the shelf. A biological ‘hotspot’ has been identified southwest of St. Lawrence 
Island between the 40- and 70-m isobaths (Grebmeier and Cooper 1995). Model results 
show that a relatively strong source of on-shelf transport occurs south of this ‘hotspot’ in 
the Zhemchug Canyon region. Upwelling of nutrient-rich deep Bering Sea water due to 
eddies and shelf-basin exchange in general could enhance biological productivity in this 
region. For example, the eddy mentioned above in Section G of Chapter V, had a mean 
salinity of 33.12, which can be converted to a nitrate concentration of 19.06 pM and also 
a silica concentration of 41.98 pM (based on relations from data in Cooper et al. 1997). 
If the 3-month (October-December 1993) mean volume transport (0.416 Sv) across 
Zhemchug Canyon section is used, this yields a total of 3,912 kg of nitrate and 8,350 kg 
of silica advected onto the Bering Sea shelf through Zhemchug Canyon due to the eddy. 
Focusing on upwelled water (deeper than 100 m), the mean salinity is 33.36, which 
would be associated with a nitrate concentration of 21.13 pM and a silica concentration 
of 45.61 pM. Over the 3-month lifetime of the eddy, the mean volume transport (for 
water deeper than 00m) is 0.11 Sv, which would result in 1,146 kg of nitrate and 2,426 kg 
silica advected onto the Bering Sea shelf. It is important to note that these relationships 
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between salinity and nutrient concentrations are rough estimates. However, to my 
knowledge this is the first calculation of the amount of nitrate or silica provided to the 
shelf via an eddy along the Bering Sea slope. 

Some of the highest silicate concentrations in the Bering Sea have been measured 
in the Gulf of Anadyr (on the western part of the shelf; Cooper et al. 1997); however 
these are diluted and biologically utilized within the euphotic zone downstream, as this 
water moves northward toward Bering Strait (Clement et al. 2005). Modeled EKE fields 
show that the northern Bering Sea maintains year-round high energy and mixing, 
especially in Bering and Anadyr straits. Notably, these regions of high EKE are found 
just upstream of highly productive areas in the Bering Sea (e.g., the Chirikov Basin and 
the region just north of Bering Strait) that have been identified in previous studies such as 
Grebmeier et al. (1988), Springer and McRoy (1993), and Grebmeier and Dunton, 
(2000). This suggests that high-nutrient Anadyr Water is mixed into the euphotic zone as 
it flows generally northward and, upon encountering a region of lower EKE, can support 
water column primary production and the settling of organic matter to the benthos. 

C. DOWNSTREAM EFFECTS 

Satellite records of the Arctic sea ice cover show a negative trend in ice extent 
and concentration since 1979 (Cavalieri et al. 1997; Comiso 2002; Meier et al. 2005; 
Parkinson et al. 1999; Serreze et al. 2003; Stroeve et al. 2005). This trend, superimposed 
over large seasonal and interannual variability (Comiso et al. 2003), has been coincident 
with the increase of the Northern Hemisphere Annular Mode (NAM also known as the 
AO) represented by a reduced winter weather regime over mid- to high-latitude 
continental regions of the Northern Hemisphere between late 1960s and mid 1990s 
(Thompson and Wallace 2001; Rigor et al. 2002). The warming as represented by the 
declining ice pack in the Arctic Ocean in summer has intensified during the late 1990's 
and into the 2000’s, especially in the western Arctic Ocean as observed from satellites 
and in situ measurements (Comiso et al. 2003; McPhee et al. 2003; Morison et al. 2002). 
The primary result of this warming has been melting of the Arctic ice pack and the 

dramatic reduction of summer sea ice cover (Maslowski et al. 2007; Rothrock et al. 1999, 
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2003; Serreze et al. 2003; Wensnahan et al. 2007). Updated estimates of Arctic 
September sea ice extent indicate a negative trend of over 10% per decade since 1979 
(Comiso et al. 2006; Stroeve and Maslowski 2007; Stroeve et al. 2005), with some of the 
lowest extents occurring in 2002-2007. During these summers, the reductions in sea ice 
cover were particularly pronounced in the western Arctic and in particular off the shores 
of Alaska and Siberia. It is important to note that this accelerated warming in the 2000s 
has occurred under a relatively neutral AO regime (while wanning has been typically 
associated with high positive AO index), which poses important questions about the 
actual role of AO in sea ice variability (Overland and Wang 2005). 

Recent reduction of the Arctic ice pack has been primarily associated with 
anomalies of surface air temperature and atmospheric circulation over the Arctic (Francis 
et al. 2005; Francis and Hunter 2006) and those in turn have been li nk ed to the AO (Rigor 
et al. 2002). Such studies typically assume that external atmospheric forcing plays the 
dominant role and neglect effects of processes internal to the Arctic Ocean. Two 
processes that tend to be especially overlooked are the oceanic thermodynamic control of 
sea ice via under-ice ablation and lateral melt at marginal ice zones and downstream 
along the advection path of warm Pacific Water (Shimada et al. 2006). The absolute 
magnitude and long-tenn variability of oceanic heat flux into the western Arctic are not 
well known from observations and are typically poorly represented in models. However, 
those ice-ocean interactions at seasonal to interannual time scales may act to decorrelate 
AO forcing, which could help explain some of the timing issues between 
AO/atmospheric forcing and sea ice variability, especially during the last decade 
(Overland and Wang 2005). 

It is hypothesized that the oceanic heat flux into the western Arctic Ocean is one 
of the main driving forces acting to reduce the sea ice cover there. In addition to the 
atmospheric forcing at the surface, the oceanic heat distributed above the halocline 
modulates the state of sea ice cover and determines regions of net growth/melt of sea ice 
and variability in multi-year and first-year ice distribution. The removal of sea ice along 
the Alaskan coast and from the Chukchi shelf for prolonged periods of time acts to 
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increase oceanic heat content in the upper ocean and its effects on sea ice in the Beaufort 
and East Siberian seas by increasing the surface albedo. Therefore, it is possible for the 
flow through Bering Strait to have an important effect of the initial melting of sea ice 
within the Chukchi Sea. 

Heat fluxes through Bering Strait and through the Chukchi shelf appear to 
influence the distribution and thickness of sea ice (Shimada et al. 2006; Woodgate et al. 
2010). Woodgate et al. (2010) combined near-bottom measurements of the velocity and 
temperature with satellite-derived SST to produce time series of the heat flux through 
Bering Strait. Estimates show a doubling of the heat flux between 2001 and 2007 (2- 
3xl0 20 J/yr to 5-6x10 20 J/yr). Woodgate et al. (2010) suggest that heat flux variability 
across Bering Strait is slightly larger than the variability of the annual shortwave input to 
the Chukchi Sea between 1998 and 2007 (~3 and 4.5x10 20 J/yr, respectively). Shimada et 
al. (2006) note that the distribution of warm Pacific Summer Water is similar to the 
spatial pattern of recent sea ice reduction since the late 1990s. However, there does not 
appear to be a direct link to Bering Strait throughflow, as increases in Pacific Summer 
Water temperatures in the basin do not correlate with the temperature of the upstream 
source water from the Bering Sea. 

Model results suggest that moderate increases in the Bering Strait heat flux in the 
early 2000s combined with solar heating on the shelf of the Chukchi Sea produced an 
increase in the northward oceanic heat flux across the Chukchi Sea. Two-thirds of the 
heat flux entering the basin was contained in the Alaska Coastal Current, which showed 
much warmer temperatures and higher velocities during September of the 2000s as 
compared to the previous two decades (Maslowski and Clement Kinney 2011). 

Freshwater flux from the Bering Sea into the Chukchi Sea is an important factor 
affecting stratification and the maintenance of the Arctic Ocean halocline (e.g., Aagaard 
et al. 1985a). This halocline separates the upper mixed layer from the deeper, warmer 
Atlantic layer. Any disturbance or reduction in the strength of the halocline could allow 
warm Atlantic Water to be mixed into the surface layer and possibly come into contact 
with the overlying sea ice. Multiple model results and observations presented in Chapter 
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VII show no long-term trends in the freshwater flux time series across Bering Strait over 
the 1979-2004 time period. However, the available measurements and multiple models 
all have problems representing the ACC contribution to the total freshwater flux (see 
additional comments below). 
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IX. OVERALL CONCLUSIONS 


The use of a high-resolution (~9km) ice-ocean model with a large domain (pan- 
Arctic) allowed for a comprehensive analysis of the circulation and property exchange 
within the Bering Sea, as well as the communication from the south and to the north. 
This model represents a large step forward in the ability to simulate the mesoscale 
activity in the North Pacific and deep Bering Sea basin, which was shown to be a strong 
control on the flow into and out of the western Aleutian Island passes and straits. A long 
integration time (1979-2004) reveals the large interannual variability in the flow through 
both the southern and northern straits of the Bering Sea. Results show that representation 
of the mesoscale features within the Alaskan Stream, Western Subarctic Gyre, and the 
deep Bering Sea basin is critical in order to simulate the communication between the 
North Pacific and the Bering Sea. In addition, long-tenn continuous observations of the 
flow through the western passes are necessary to better understand the actual variability 
there and to constrain models. 

The Bering Slope Current appeared to be more a system of eddies rather than a 
continuous feature, which indicates the need for a fully eddy-resolving model to represent 
the complex dynamics along the slope. There does not appear to be a north-south split of 
the BSC at Cape Navarin, instead the model results show that the main source water for 
the Gulf of Anadyr comes from the outer shelf flow, roughly between the 100 and 200-m 
isobaths. Important areas of shelf-basin exchange were identified along the Bering Slope, 
particularly within the canyons, such as Zhemchug and Bering canyons. Peaks in the on- 
shelf transport of salt and heat were shown to be associated with cyclonic eddies in 
Zhemchug Canyon, which makes this area important for supplying nutrients to the Bering 
Sea shelf, as well as the entire Western Arctic. An eddy was estimated to supply 3,912 
kg of nitrate and 8,350 kg of silica during the 3 months when it was present within the 
canyon. 

The model represented a major anomaly in the flow and sea ice distribution across 
the northern Bering Sea shelf. The anomaly, which also was apparent in satellite-derived 
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sea ice concentration data, was caused by an unusual wind pattern. Strong winds from a 
more easterly direction created an Ekman transport on the shelf which reversed the 
Anadyr Current and greatly increased the northward flow through Shpanberg Strait. 
Additional model results within the northern Bering Sea showed that increased EKE in 
Bering and Anadyr straits may help explain the areas of high biological productivity 
located just downstream in the Chirikov Basin and north of Bering Strait. 

Estimates of the transport through Bering Strait from multiple models ranged 
from 0.67 (+/- 0.03) to 1.29 (+/- 0.06) Sv in the mean, compared to the observational 
estimate of 0.8 +/-0.2 Sv. Thus, most of the models agreed with the observational 
estimate to within errors. However, there is still a need to reduce the uncertainty in the 
accurate estimation volume transport and property fluxes through Bering Strait, both 
from observations and model results. Similar to the need for more observations in the 
western Aleutian passes, there is a need to observe the horizontal and vertical structure of 
the flow through Bering Strait, including the western portion, on a continuous basis. An 
effort is currently underway to provide more observations within Bering Strait and these 
measurements will be important for many future studies of the Pacific-Arctic region. 
Increases in model resolution and improvements in model bathymetry should improve the 
representation of the flow structure in Bering Strait, as indicated by initial spin-up results 
from a new, higher resolution version of the NAME model. 
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X. FUTURE DIRECTIONS 


The lack of tides in the NAME model limits its ability to properly represent short 
period (-“1-15 days) variability in areas where tides are important, such as within the 
Aleutian Island Passes (e.g., Kowalik 1999; Stabeno et al. 2002). Residual tidal currents 
occur when the energy of the basic tidal constituents is transferred to shorter and longer 
periods and the oscillating velocity interacts either with the sea surface or with bottom 
topography. Although observations on the formation of new wave periods are very few, 
it appears that residual currents can be an important contribution to the mean flow 
(Kowalik 1999). 

Model bathymetry representation of some features at 9-km grid cell spacing, 
namely important passes and straits, does not always agree with the observed bathymetry. 
In general, there are local inaccuracies in the digital bathymetric dataset and the 
resolution of the digital bathymetry is coarse relative to the dimensions of the Aleutian 
Island passes and straits of the northern Bering Sea. Increased model resolution and an 
improved digital bathymetry for the passes and straits are expected to improve the model 
skill at these locations. Improvements in the 2-km version of the NAME model 
bathymetry should rectify some of these problems. 

It seems likely that the seasonally- and surface-intensified Alaska Coastal Current 
(ACC) volume transport is not fully accounted for in large-scale models due to spatial 
resolution limitations. At the same time, the estimates presented here from observations 
also lack continuous measurements in the surface layers and near the coast. Although 
estimates of the contributions from the ACC and stratification have been made by 
Woodgate et al. (2006, 2010), quantification of the seasonal contribution to the Bering 
Strait transport is yet to be computed from either observations or models. The freshwater 
flux, which has a significant influence on the density structure of the Arctic Ocean (e.g., 
Aagaard et al. 1985b), would also be better measured if salinity infonnation could be 
obtained in the upper layers and nearby the coast. Similarly, for heat flux, it is crucial to 
get information on the upper layers where maybe 1/3 of the heat is advected [see 
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Woodgate et al. (2010) who used satellite-derived sea surface temperatures to estimate 
the contribution from the upper layers]. An international effort is currently underway 
with 8 moorings placed in the Bering Strait region. New information from these 
moorings will be important for better understanding details of the flow through the strait 
and for constraining models. 

While it is encouraging that, in many of the larger-scale models, fluxes of volume, 
heat and salt through the Bering Strait are of the right order of magnitude and in 
interannual terms show correlated variations with observations, there are still 
discrepancies. These have to be considered when using model results to look at the role 
of Pacific waters in the Arctic. There is also a need for models with higher spatial 
resolution in the Bering Strait region. The ACC is only order 10 km in width and thus 
not resolved by global or regional Arctic models with resolutions of 7-22 km (Table 7). 
The implementation of higher-resolution (2 km or less) regional models should improve 
estimates of the volume, heat and freshwater fluxes in the strait. The challenge is to be 
able to capture small-scale features, such as the Alaska Coastal Current and mesoscale 
eddies in the strait itself and its immediate vicinity. The modeling community is working 
toward that goal to properly represent such features. 
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